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PREFACE 

This  final  report  presents  the  results  of  work  performed  by  the 
Lockheed  Missiles  & Space  Company,  Inc.,  Huntsville  Research  & Engineer- 
ing Center,  under  Contract  DOT-TSC-1098  for  the  Department  of  Trans- 
portation Systems  Center,  Cambridge,  Mass.  The  objective  of  the  contract 
was  to  modify,  calibrate,  and  deploy  a mobile  laser  Doppler  velocimeter 
system  for  the  measurement  of  wake  vortices  and  winds  at  terminal  areas. 
The  period  of  performance  for  this  study  was  from  July  1975  through  June 
1976.  Lockheed-Huntsville  personnel  contributing  to  this  effort  were  E.  W, 
Coffey,  C.  E.  Craven,  B.  B.  Edwards,  E.W.  Feese,  E.J.  Gorzynski,  J.  L. 
Jetton,  A.  J.  Jordan,  M.  C.  Krause,  G.  M.  Miller  and  K.  R.  Shrider.  The 
Contracting  Officer's  Representative  for  this  work  was  Dr.  J.N.  Hallock. 
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1.  INTRODUCTION  AND  SUMMARY 


1.1  BACKGROUND 

Considerable  effort  is  currently  being  devoted  to  the  development  of 
instrumentation  to  remotely  sense  atmospheric  flow  phenomena.  Some  of 
the  avenues  being  pursued  are  both  active  and  passive  acoustic,  optical,  and 
radio  methods.  A useful  survey  of  such  methods  is  presented  in  Ref.  1.  Two 
advantages  of  remote  sensors  are  that  flow  conditions  can  be  ascertained  in 
regions  of  space  where  it  would  not  be  convenient  to  locate  conventional  in- 
strumentation, and  no  modification  of  the  flow  at  the  point  of  interest  is  intro- 
duced by  their  use.  The  laser  Doppler  velocimeter  (LDV)  is  a particularly 
attractive  device  for  remote  sensing  of  atmospheric  phenomena.  In  the  LDV 
system,  laser  radiation  is  backscattered  from  moving  particulates  in  the 
atmosphere  and  is  used  to  establish  the  velocity  of  the  flow.  Since  optical 
tracking  of  the  laser  focal  volume  is  possible,  a scanning  LDV  system  can 
rapidly  determine  the  velocity  field  over  a large  region  in  space.  A CO^ 
laser  Doppler  velocimeter  system  possesses  the  following  advantages 
over  other  remote  sensing  techniques:  (1)  the  sensing  volume  can  be  varied 
with  ease  as  only  optic  pointing  and  focusing  operations  are  involved;  (2)  the 
ambient  aerosol  provides  an  adequate  scattering  target;  and  (3)  the  sensing 
mechanism  is  non-mechanical  which  results  in  the  potential  for  a high  fre- 
quency turbulence  sensor. 

The  feasibility  of  utilising  a LDV  system  for  the  remote  sensing  of 
boundary  layer  winds  and  for  the  detection,  tracking  and  measurement  of 
aircraft  wake  vortices  has  been  demonstrated  (Refs. 2,  3 and  4).  However, 
the  development  of  an  effective  LDV  system  for  monitoring  wind,  wind  shear, 
and  wake  vortices  at  terminal  areas  required  further  refinement  and  appli- 
cation of  this  technology  Including  the  following  tasks: 


ft 


Design  and  fabrication  of  a compact,  mobile,  self-supporting  LDV  « 

system. 

Improvements  to  the  automatic  optical  scanner,  display,  and  soft- 
ware and  accommodation  of  both  wind  and  vortex  tracking  modes 
of  operation. 

Comprehensive  field  testing  of  the  LDV  system  in  both  wind  and 
vortex  tracking  modes  to  establish  the  baf  c operational  capabilities, 
resolution,  and  integrity  of  the  system. 

This  technical  report  deals  with  the  above  tasks  and  summarizes  the  efforts 
carried  out  by  Lockheed -Huntsville  to  develop  a wind,  wind  shear,  and 
wake  vortex  remote  sensor.  Tasks  1 and  2 were  undertaken  by  Lockheed- 
Huntsville  using  Company  funds.  A mobile  self-supporting  LDV  system,  the 
Lockheed -Huntsville  LDV  van,  with  automatic  scanning  capabilities  for  wind 
and  wake  vortex  measurements  was  designed,  fabricated,  and  field  tested  in 
1975.  Based  on  the  successful  demonstration  of  the  LDV  system,  the  Depart- 
ment of  Transportation,  Transportation  Systems  Center,  contracted  Lockheed- 
Huntsville  to  make  further  refinements  to  the  mobile  LDV  system  and  to  deploy 
it  at  the  John  F.  Kennedy  International  Airport  (JFK)  and  to  evaluate  its  opera-  ' 

tional  capabilities.  The  results  of  the  research  and  development  effort  are 
the  subject  of  this  technical  report. 

1.2  PROGRAM  OBJECTIVES 

The  research  program  focused  on  the  evaluation  of  the  Lockheed- 
Huntsville  LDV  for  providing  wind,  wind  shear,  and  wake  vortex  measure- 
ments in  support  of  airport  operations.  The  program  encompassed  the 
following  tasks; 

Completion  of  minor  modifications  to  the  Lockheed -Huntsville  LDV 

to  enable  tracking  of  vortices  over  a long  period  of  time  : 

Calibration  of  the  laser  system  in  its  various  scamning  modes 
Installation  of  displays  to  monitor  the  operation  of  the  system  on-line 

Collection  of  a data  base  of  both  wind  and  wake  vortex  measurements 

at  JFK  f 

Establishment  of  the  overall  performance  capabilities  of  the  system 
for  wind,  wind  shear,  and  wake  vortex  measurements  based  on  the 
analysis  of  the  above  measurements. 
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To  achieve  these  program  objectives,  the  LDV  system  was  deployed  at  two 
test  sites:  at  the  Huntsville  Jetplex  in  Huntsville,  Alabama,  and  at  the  John 
F.  Kennedy  (JFK)  International  Airport  in  New  York.  The  purpose  of  the 
Huntsville  Jetplex  tests  was  to  calibrate  the  system  and  to  obtain  sample 
wake  vortex  measurements.  Following  the  successful  demonstration  of  the 
system  at  the  Jetplex,  the  LDV  system  was  transported  to  JFK.  During  the 
JFK  field  tests,  spanning  30  days,  the  LDV  observed  and  recorded  wake 
vortex  trajectories  during  normal  landing  operations.  In  addition,  wind  field 
measurements  were  carried  out  with  the  LDV  and  compared  with  JFK  meteor- 
ological tower  measurements.  The  report  discusses  the  check  out,  calibra- 
tion, and  operation  of  the  mobile  laser  Doppler  system  during  the  two  tests. 


1.3  REPORT  FORMAT 

A discussion  of  the  mobile  laser  Doppler  system  is  presented  in  the 
following  sections.  The  development  of  the  LDV  system  is  addressed  first 
in  Section  2 followed  by  a description  of  the  computer  software  algorithms 
in  Section  3.  The  field  tests  of  the  LDV  are  discussed  in  Section  4 and  the 
wake  vortex  measurements  are  presented  in  Section  5.  Wind  and  wind  shear 
measurements  are  discussed  in  Section  6.  A summary  and  overview  of  the 
program  is  given  in  Section  7. 


2.  LASER  DOPPLER  SYSTEM  DEVELOPMENT 

2.1  SYSTEM  DESCRIPTION 

A remote  sensing  system  has  been  developed  by  Lockheed-Huntsville 
for  tracking  wake  vortices  and  for  measuring  a three-dimensional 
wind  profile.  The  basic  hardware  consists  of  a laser  and  associated  optical 
systems,  the  scanning  system,  and  the  display/processing  system  as  shown 
in  Fig.  2-1.  The  apparatus  is  housed  in  a mobile  van,  sketched  in  the  upper 
right-hand  corner  of  Fig.  2-1.  A description  of  the  LDV  system  includes 
discussion  of  the  laser  Doppler  principle,  optical  system,  scanning  system, 
signal  processing  system,  and  data  recording  and  display  system. 

2.1.1  Laser  Doppler  Principle 

An  LDV  wind/vortex  sensor  involves  measurement  of  the  Doppler  spec- 
trum of  laser  radiation  backscattered  by  atmospheric  aerosols.  The  instru- 
ment must  incorporate  means  to  transmit  the  laser  radiation  to  the  region  of 
interest,  collect  the  radiation  backscattered  from  the  atmospheric  aerosol, 
and  to  photomix  on  a photodetector  the  scattered  radiation  and  a portion  of 
the  traomitted  beam.  A variable  frequency  component,  at  the  Doppler  shift 
frequency,  is  generated  at  the  photodetector  which  is  translatable  into  an 
along -optic -axis  wind  velocity  component  using  appropriate  electronics.  The 
magnitude  of  the  Doppler  shift,  Af,  is  given  by  the  equation 

Af  = I lv|  cose, 

where 

V s velocity  vector  in  the  region  being  sensed 

A s the  laser  radiation  wavelength  (10.6  microns  for  the 
CO2  laser),  and 

e - the  angle  subtended  by  the  velocity  vector  and  the 
optic  system  line  of  sight. 


2-1 


SCANNING  SYSTEM  MINICOMPUTER  AND  DISPLAY 


Using  the  CO^  laser,  1 m/sec  of  line-of-sight  velocity  component  yields  a 
Doppler  shift  of  188  kHz.  Thus,  measurement  of  the  Doppler  shift  frequency, 
Af,  yields  directly  the  line-of-sight  velocity  component  | V|  cosO.  Some 
advantages  of  the  laser  Doppler  method  are;  (1)  the  Doppler  shift  is  a di- 
rect absolute  measure  of  the  velocity  (for  example,  a hot  wire  yields  velocity 
via  a cooling  effect  on  the  wire),  (2)  the  ease  with  which  the  sensing  volume 
can  be  varied  (optics  pointing  and  focusing  operations  only  being  involved); 

(3)  the  ambient  aerosols  provide  sufficient  scattering,  thus  enabling  opera- 
tion in  clear  air  conditions;  and  (4)  the  ambient  aerosol  tracer  has  a small 
inertia  and  responds  quickly  to  variations  in  windspeed  and  could  thus  be  a 
good  turbulence  indicator. 

A useful  instrument  must  incorporate  means  to  scan  the  system's  reso- 
lution volume  in  a prescribed  manner  and  to  effect  the  required  signal  process- 
ing, on-line  read-out,  and  permanent  recording  requirements.  The  hardware 
implementation  of  the  field  laser  Doppler  unit  used  during  the  investigation 
is  discussed  in  the  following  subsections. 

2.1.2  Optical  System 

The  optical  system  is  of  a monostatic  design  and  uses  a continuous 
wave  laser.  The  system  depends  on  focusing  the  transmitter  telescope  at  the 
location  of  interest  for  its  spatial  resolution  property.  Details  of  the  optical 
arrangement  are  shown  in  Fig.  2-2. 

A horizontally-polarized,  15-watt,  continuous  wave  CO2  laser  beam 
(10.6  micron  wavelength)  emerges  from  the  laser  (^and  is  deflected  90 
degrees  first  by  a mirror  ^^and  then  by  a 90%  reflecting  beamsplitter 
The  approximately  0.23  in.  diameter  beam  then  passes  through  a Brewster 
window  ^^and  a CdS  quarter  waveplate  which  converts  the  beam  to  cir- 
cular polarization.  The  beam  impinges  on  the  secondary  mirror  C^and  is 
expanded  and  reflected  into  the  primary  mirror  (1  ft  diameter)  and 
then  focused  out  into  the  atmosphere.  A wire  stop  eliminates  most  of 
the  secondary  mirror  reflection  of  the  outgoing  beam.  A small  portion  of 
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the  original  laser  beam  is  transmitted  through  the  beamsplitter  (^and  is 
used  as  a local  oscillator  after  being  rotated  to  vertical  polarization  byr  a 
half  waveplate  Energ/  scattered  b^  aerosols,  at  the  focal  volume 

is  collected  by  the  primary  mirror  collimated  by  the  secondary  and 
passed  through  the  quarter  waveplate  ^).  The  quarter  waveplate  changes 
the  polarization  of  the  aerosol  backscattered  radiation  from  circular  to  vertical 
linear  polarization.  The  vertically  polarized  beam  is  approximately  78%  re- 
flected off  the  Brewster  window  and  directed  via  a mirror  through  the 
beamsplitter  where  it  is  combined  with  the  local  oscillator  radiation. 

After  passing  through  the  collecting  lens  the  two  beams  are  photomixed 
on  the  detector  in  a heterodyne  configuration.  The  electrical  output  of 
the  detector  is  amplified  with  a 5 MHz  bandwidth,  20  dB  gain  low  noise 
type  preamplifier  and  fed  into  a spectrum  analyzer 

An  alternative  operating  configuration  consists  of  using  the  portion 
of  the  outgoing  beam  backscattered  into  the  interferometer  by  the  secondary 
mirror  as  the  local  oscillator  beam.  This  mode  of  operation  is  less 
susceptible  to  optical  misalignment  difficulties  and  was  the  technique  used 
during  most  of  this  investigation.  When  incorporated,  the  optical  leg 
was  deactivated  by  the  beam  stop  (^and  the  wire  stop  removed. 

2.1.3  Optical  Scanning  System 

To  give  the  flexibility  needed  to  operate  in  various  modes  of  operation 
(Section  2.3)  the  anteana  arrangement  shown  in  Fig. 2- 3 is  used.  The  mirror 
arrangement  AB  can  be  rotated  about  the  vertical  axis,  thus  producing  the 
Velocity -Azimuth  Display  (VAD)  or  conical  scan  mode  of  operation.  Mirror 
A is  adjusted  to  deflect  the  beam  into  a plane  normal  to  the  plane  of  the  paper 
and  at  a zenith  angle  corresponding  to  the  required  conical  scan  angle. 

A plMs  or  an  arc  in  space  can  be  interrogated  by  using  the  system's 
elevation  scan  capability  which  consists  of  rotating  mirror  A about  a horizontal 
axis.  The  transfer  of  this  capability  from  mirror  C to  mirror  A and  the  in- 
crease of  the  elevation  scan  angle  range  from  60  to  180  degrees  was  one  of 


Conical  Scan 
Rotation  Axis 
(Vertical) 


I Fig.  2-3  - Scan  Equipment  Added  to  LDV 
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the  tasks  of  this  program.  This  added  capability  admitted  considerable 
flexibility  ip  the  tracking  of  wake  vortices.  The  scanning  hardware  is  de- 
ployed on  the  mobile  van  as  shown  in  Fig.  2-4. 

The  system's  focal  volume  is  range  scanned  by  varying  the  separation 
of  the  telescope  secondary  mirror,  E,  and  the  primary  mirror,  D (Fig.  2-3). 
This  is  done  by  varying  the  position  of  the  mirror,  E,  in  a controlled  manner 
by  an  electric  motor/optical  encoder  combination.  The  scanning  system  is 
addressed  by  a control  panel  incorporating  thumbswitch  controls  and  LED 
monitors.  The  system's  scan  capabilities  are  summarized  in  Table  2-1  and 
Fig. 2-5. 


2.1.4  Signal  Processing  System 

The  LDV  signal  processing  system  consists  of  a detector,  spectrum 
analyser,  frequency  voltage  converter  (tracker),  input/output  coupler,  and 
peripheral  hardware.  A block  diagram  of  the  basic  signal  processing  system 
is  shown  in  Fig. 2-6.  The  heterodyned  laser  rettirn  is  converted  into  an  elec- 
trical signal  by  the  detector  and  is  amplified.  The  frequency  content  of  the 
signal  is  determined  by  the  spectrum  analyzer.  The  resulting  signal  is  dis- 
played and  its  intensity  and  frequency  ultimately  recorded  in  digital  form. 


The  Doppler  content  of  the  photodetector  output  is  extracted  through  the 
use  of  a sampled  spectrum  analyzer  which  provides  frequency  spectra  at  a 
70  Hz  rate.  A typical  Doppler  wind  spectrum  is  shown  in  Fig.  2-7.  To  yield 
a line-of- sight  velocity  estimate,  a voltage  is  generated  which  has  the  same 
time  behavior  as  the  Doppler  shift,  f^,  as  given  by  the  peak  of  the  spectrum. 

The  implementation  of  this  technique  is,  in  essence,  a recursive  com- 
parison method.  The  spectrum  analyzer  scan  is  driven  by  a sawtooth  voltage 
derived  from  a D/A  converter,  the  input  to  which  is  counter  clocked  at  a con- 
stant rate,  hence  the  digital  number  output  of  the  counter  represents  frequency 
on  a linear  scale.  At  each  new  count,  the  spectrum  analyzer  output  is  con- 
verted to  a digital  representation  by  an  A/D  converter,  and  the  binary  number 
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Table  2-1 

SCAN  CAPABILITY 


RANGE 

Maximum  Limit:  2600  ft 

Minimum  Limit:  52  ft 

Scan  Frequency:  0.1  to  6.9  Hz 

ELEVATION 

Maximum  Limit;  90  deg 

Minimum  Limit:  0 deg 

Hard  Limit;  3 deg  with  Override 

Scan  Frequency;  0.1  to  0.5  Hz 

VAD  MODE 

Measurement  Altitude:  33  to  2100  ft 

Measurement  Time/ Altitude;  5 sec 

Sample  Rate;  1 to  7 Cycles 

Number  of  Altitudes;  8 

MULTIMODE 

Elevation  Coverage:  3 to  90  deg  Upwind  and 
Downwind 

Scan  Plane  Azimuth:  360  deg 

Vertical  Line  Scan;  52  to  2100  ft 

Overhead  Arc  Scan;  90  deg  Coverage  Maximum 

ACCURACY 

Range:  1 ft  at  98  ft,  98  ft  at  984  ft 

Elevation:  0.25  deg 

es 


Fig.  £.-b  ~ LDV  Signal  Processing  System 


0.75  MHz  2 MHz 

(~  13  ft/sec)  (34.8  ft/sec) 


Frequency  (velocity) — ^ 


Fig.  2- 7 - Typical  LDV  Wind  Signature  as  Displayed  by  Spectrum  Analyzer 
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representing  the  current  sample  is  compared  with  that  obtained  on  the  pre- 
vious count.  If  the  current  one  is  the  larger  of  the  two,  it  is  saved  by  storing 
in  a latch,  along  with  the  binary  number  representing  its  frequency;  if  it  is 
smaller,  the  previous  one  is  retained  until  the  next  comparison.  This  process 
is  continued  for  the  entire  sweep.  It  is  evident  that  the  number  remaining  in 
the  frequency  store  latch,  when  the  sweep  is  completed,  corresponds  to  the 
highest  signal  power  observed,  i.e.,  the  peak  of  the  spectrum.  At  the  end  of 
each  sweep,  the  new  peak  frequency  found  replaces  that  obtained  on  the  pre- 
vious sweep.  An  example  of  the  output  is  shown  in  Fig.  2-8  for  tl  e case  of  an 
FM  signal  of  center  frequency  2.0  MHz  (f^)  modulated  to  +200  kHz  about  f^  at 
a 5 Hz  rate  sinusoidally. 


f =5  Hz  Sine  Wave 
m 


Oscilloscope  Data 

Horiz.=  O.l  sec/div 
Vert.  = 1 V/div 


Fig.  2-8  - Output  of  Signal  Processor  for  FM  Modulated  Input 


A provision  is  included  for  tracking  single  sideband-suppressed  carrier 
signals,  with  an  identification  of  upper  or  lower  sideband  such  that  if  used 
in  conjunction  with  an  acousto-optic  modulator  the  unit  can  discriminate  the 
sign  of  the  Doppler  shift.  The  signal  feedthrough  at  the  translated  frequency 
can  also  be  discriminated  against  digitally,  thus  elin.inating  the  need  for  a 
"notch  filter." 
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The  raw  spectral  information  (output  of  the  spectrum  analyzer)  is  also 
made  available  to  the  Systems  Engineering  Laboratories  (SEL)  810A  data 
logging  minicomputer  which  is  programmed  to  generate  its  own  estimate  of 
the  spectral  peak  and  is  discussed  later  in  Section  3. 

2.1.5  Data  Recording  and  Display 

The  primary  data  gathering  function  is  performed  by  a SEL  810A  general 
purpose  minicomputer.  Data  acquired  by  the  Mobile  Atmospheric  Unit  is  for- 
matted by  the  computer  software  and  stored  on  magnetic  tape  for  subsequent 
processing.  The  SEL  7-track  tape  control  and  magnetic  tape  units  allow 
digital  recording  of  data  at  800  bpi  at  45  ips.  The  data  logged  by  the  com- 
puter includes: 

a.  All  scan  volume  location  parameters 

b.  "Mode  of  operation"  identifier 

c.  The  instantaneous  line-of-sight  velocity  information 

d.  The  Doppler  spectrum  peak  strength 

e.  Full  spectrum  intensity  and  frequency  information 
(optical) 

f.  A data  quality  identifier. 

Properties  of  the  Doppler  spectrum,  namely  the  amplitude  and  frequency 
corresponding  to  the  sp>ectral  peak,  are  obtained  as  a result  of  on-line  com- 
puter processing;  the  frequency  is  also  obtained  by  the  spectral  peak  locator 
(velocity  processor)  discussed  previously.  The  latter  allows  some  flexibility 
for  on-line  operator  displays  (see  below). 

The  velocity  processor  estimate  of  the  instantaneous  line-of-sight 
velocity,  updated  at  a 70-Hz  rate,  is  available  in  analog  format  which  can 
be  recorded  directly  on  a strip  chart  recorder  (an  option  which  is  extremely 
useful  during  the  VAD  mode  of  operation  for  monitoring  the  characteristic 
profile).  During  vortex  tracking  operations  the  velocity  processor's  output 
j is  used  to  modulate  the  intensity  of  a CRT  beam  which  is  driven  to  trace 

I on  the  tube  face  an  exact  replica  of  the  scan  pattern  generated  by  the 
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scanning  system.  The  display  shows  on-line  the  location  of  high  velocity 
regions.  For  vortex  tracking  purposes,  the  display  aids  the  operator  in 
selecting  the  optimum  scan  adjustments  so  that  the  trajectory  of  the  wake 
vortex  can  be  recorded.  The  integration  of  such  a vortex  display  into  the 
system  was  a required  task  in  the  Scope  of  Work.  The  information  in  this 
subsection  is  supplemented  by  the  discussion  in  Section  2.4  relating  to  data 
transfer  and  formatting. 

Some  overall  views  of  the  mobile  unit  hardware  as  used  during  the 
program  are  shown  in  Figs.  2-9  through  2-12. 

2.2  SYSTEM  SPECIFICATIONS 

Specifications  for  the  subsystems  and  components  are  as  shown  in 
Table  2-2. 

2.3  MODES  OF  SYSTEM  OPERATION 
2.3.1  Winds  Aloft  Sensing 

Using  the  basic  system  outlined  previously  it  is  possible,  by  scanning 
operations,  to  determine  the  three-component  wind  field  at  any  specified 
altitude  between  50  and  2000  feet.  Since  the  LDV  system  measures  the  wind 
velocity  at  a given  point  in  space  along  the  optical  system's  line  of  sight,  it 
is  necessary  to  sample  the  line- of- sight  velocity  at  different  points  in  space 
to  compute  the  full  three-dimensional  wind  field.  The  scanning  method  em- 
ployed  is  commonly  referred  to  as  the  Velocity  Azimuth  Display  (VAD)  tech- 
nique which  was  used  by  Lhermitte  and  Atlas  in  conjunction  with  a micro- 
wave  radar  (Ref.  5). 

The  telescope  is  focused  at  the  altitude  of  interest,  the  beam  being  di- 
rected at  a zenith  angle,  a.  The  beam  is  then  scanned  in  azimuth,  thus  tracing 


Also  known  as  conical  scan  technique  because  of  beam  scanning  configuration. 
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Fig.  2-9  - LDV  Monitoring  Wake  Vortex  Generated  by  L-1011  at  Huntsville  Airport 


rst  Rack,  Computer  in  Second  Rack,  Digital  Tape  Unit  Aft,  and  Optics 
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Table  2-2 

LOCKHEED  DETAILED  COMPONENT  SPECIFICATIONS 


Overall  Performance  Specifications 

Range  « 32  ft  to  2600  ft  Line  of  Sight 

Velocity  Measurement  Threshold:  ^ -7  ft/sec  Line  of  Sight 

Measurement  Accuracy:  Velocity  » + 1.64  ft/sec,  Direction  « ^3  deg. 

1. 

Laser 

Type:  Honeywell  Model  9000  with  Recharging  Capability 

Power:  15  W at  P(20)  (Nominal) 

Mode:  TEMgcn*  (Single  Mode) 

Stability:  Long  Term,  +2%;  Short  Term,  ;f0.5% 

Beam  Diameter:  0.019'7~ft 

2. 

Detector 

Type:  Rockwell,  Pb-Sn-Te  P*hotovoltaic 

Quantum  Efficiency  >■  40% 

Freq.  Response:  r>  SS  MHz 

Element  Size:  x 10“'  ft^ 

Dynamic  Impedance:  >200 

Number  of  Elements:  4 

3. 

Telescope 

Diameter:  12  in.  Primary,  0.5  in.  Secondary 

Spatial  Resolution:  4 125  ft  at  1000  ft  range  at  50%  power  point 

4. 

Interferometer 

Type:  Mach-Zehnder 

Polarization  Input:  Horizontal 

Polarization  Output:  Circular 

Polarization  Recombined:  Vertical 

Components:  Mirrors  (4) 

Half  Wave  Plate  (1) 

Quarter  Wave  Plate  (1) 

Brewster  Window  (I) 

Beamsplitter  (2) 

Lens  (1) 

Aperture  (1) 

.A 


2.20 


T^ble  2-2  (Continued) 


Range/ Altitude  Scan 

Excursion:  52  ft  to  2600  ft 

Rate  (Auto):  0.1  to  6.9  Hz  in  0.1  Hz  Steps 

Rate  (Manual):  Slew  — 0.2  Hz 

Step—  6.56  ft  Increments 
Accuracy:  Static  — 0.7%  at  1148  ft  Range 

Dynamic  — 0.5%  at  561  ft  Range 

Modes:  Manual  — Slew 

Manual  — Single  Increment 
Auto  — Range 
Auto  — Altitude 

Altitude  Limits:  33  ft  to  2100  ft 

Altitude  Steps:  8 maximum,  Any  Combination 

Incremental  Altitudes  for  8 Equal  Steps 
Readout  (LED):  Range,  Altitude,  Azimuth  Angle,  Altitude  Increment 


6.  VAD  Scan 

I Measurement  Time/Alt.:  5 sec 

I Frame  Rate:  40  sec 

Cone  Angle:  0 to  60  deg.  Manually  Adjustable 


Processing  System 

Signal  Processor:  HP  Spectrum  Analyzer,  Velocity  Processor  (LMSC) 

Sample  Rate:  70  Hz 

Modes:  Translate  or  Non-Translate 

Velocity  Peak  or  Velocity  Max, 

AGC:  Auto  Threshold 
Output:  0 to  5 Vdc 
Resolution:  0.328  ft/sec 


8.  Minicomputer 

SEL  810  Data  Logging  System 
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Table  Z~2  (Concluded) 


9. 

Van  Support  Equipment 

15  kVA  4 Cycle  Generator 
(4.3  X 1.6  X 2.3  ft,  800  lb) 

24,000  Btu  Air  Conditioning 

5 kVA  Heating  (Electric) 

Leveling  Jacks 

10. 

Optional  Scan 

Elevation  Scan  Added  to  Provide  "Finger-Scan"  Mode  for  Vortex 
Tracking. 

Excursion:  0 to  90  deg 

Rate  (Auto):  0.1  to  0.5  Hz  in  0.1  Hz  Steps 

Rate  (Manual):  Slew  — 0.2  Hz 

Step—  0.5  deg  Increments 

Accuracy;  Static  — 0.8%  at  60  deg  Angle 

Dynamic  — 1.6%  at  90  deg  Angle 

Modes:  Manual  — Slew 

Manual  — Single  Increment 

Auto  — Elev. 

Readout  (LED):  Elevation  Angle 

Van  Support  — Equipped  with  air  conditioning,  heat,  power 

Equipment  distribution,  100  ft  cables,  lights,  leveling 

jacks  and  15  kVA  auxiliary  power  unit 

Signal  Processing  — Spectrum  Analyzer,  Velocity  Processor 

Data  Processor  — SEL  810A  Data  Logger 

Scanner  —Auto  Range,  Altitude,  Azimuth,  for  VAD, 

etc;  Elevation  Scan  Optional 

Operational  Equipment  — Time  Code  Generator;  Oscilloscope;  Power 

Meter;  Tesla  Coil;  and  Calibration  Target 
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out  a circle  at  the  selected  altitude  (Fig. 2-13).  The  instantaneous  mean 
radial  velocity  within  the  sensing  volume  as  measured  by  the  radar,  v^,  is 
given  by 

V = V.  sina  co8(0  - 0 ) + w cosa  , 
r h o 


Vj^  and  6^,  respectively,  being  the  speed  and  direction  of  the  horizontal  wind 
motion  and  w the  vertical  motion  at  the  height  being  sampled.  From  the 
amplitude,  phase  and  dc  component  of  the  VAD  signature  it  is  possible 
to  compute  the  horizontal  speed  and  direction  and  vertical  velocity  of  the 
wind  field. 

In  the  present  mode  of  operation  it  is  the  absolute  value  of  v^ 
is  sensed,  a deficiency  that  allows  an  ambiguity  of  r in  the  wind  direction  esti- 
mate. The  modifications  that  this  makes  to  the  VAD  profile  and  the  respective 
equations  for  the  velocity  components  are  indicated  in  Fig.  2-14. 

While  operating  in  the  VAD  mode  the  system  is  capable  of  interrogating  n 
(n  = 1 through  8)  altitudes  (that  can  be  dialed  in  by  using  thumb  switches)  in 
sequence  over  a total  time  period  of  5np  sec  where 


5 sec  = time  for  conical  sweep  at  one  altitude 

n = number  of  altitudes  to  be  interrogated 

p = number  of  VAD  scans  at  each  altitude 
(can  be  chosen  to  be  1 through  7). 

During  this  investigation  n = 8,  p = 1 were  utilized  thus  allowing  the  interro- 
gation of  eight  altitudes  every  40  seconds  and  the  sequen'^e  was  repeated  for 
the  required  test  duration. 


2.3.2  Coordinated  Range  and  Elevation  Scanning  and  Point  Sensing  Mode 
of  Operation 


The  basic  LDV  has  the  capability  of  determining  the  wind  velocity  com- 
ponent along  the  optical  system's  line  of  sight  and  within  the  system's  sensing 
volume  (Section  2.5).  The  hardware  is  arranged  to  range  scan  along  a selectable 


1 

* 

A 
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ig.2-13  - VAD  Principle  of  Operation 


Vh  = (Vj  + Sina 


Fig.  2-14  - Azimuth  Angle  Dependence  of  Measured  Velocity  Component 

line  of  sight  by  varying  the  position  of  the  telescope  secondary  mirror,  E, 
relative  to  the  primary,  D,  (Fig.  2-3)  so  as  to  scan  the  system's  sensing  volume 
linearly  between  fixed  range  limits.  The  relationship  between  the  separation, 
DE,  and  the  range  to  the  sensing  volume  is  given  by  the  lens  equation. 

This  procedure  will  result  in  the  line-of-sight  velocity  profile  along  the  chosen 
line  of  sight  and  between  the  range  limits.  The  elevation  angle  can  also  be 
scanned  in  a controlled  manner  by  scanning  the  output  mirror  shown  in  Fig.  2-3. 
A coordinated  range  and  elevation  scan  was  the  basic  mode  of  operation  during 
vortex  tracking.  Fast  range  and  slow  elevation  scans  generate  a "finger  scan" 
pattern  in  the  plane  of  interest.  Any  location  of  interest  could  be  "dialed  in" 
using  the  scanner  controls  and  that  region  interrogated  continuously.  The 
scan  patterns  traced  out  in  the  two  modes  of  operation  are  shown  in  Fig.  2-15. 

2.4  ON-LINE  DATA  MANIPULATION 

The  data  from  the  LDV  system  are  transferred  to  the  computer  along  four 
digital  input  channels.  It  is  reformatted  and  recorded  on  7-track  magnetic 
tape  for  off-line  data  processing.  Software  to  locate  the  peak  intensity  within 
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Fig.  2-15  - Scan  Configurations 


each  spectrum,  subject  to  the  constraints  of  frequency  and  amplitude  thres- 
holds and  the  associated  sample  count  number  (velocity),  is  also  incorporated. 
Additional  details  of  the  latter  process  are  included  in  Appendix  B. 

Recorded  data  consist  of  four  16-bit  words  for  each  spectrum  analyzer 
sweep.  These  four  words  correspond  to  the  four  SEL  logical  input  channels 
(60,  61,  62  and  63)  read  for  the  spectrum  analyzer  sample  count  at  which  the 
peak  intensity  occurred  during  a given  sweep.  The  logger  identifies  the  peak 
intensity  condition  independent  of  the  "data  acceptable"  condition;  four  16-bit 
words  (identifying  sample  count  location,  range  or  altitude,  angle,  and  appro- 
priate scanning  parameters)  are  saved  for  each  spectrum  analyzer  sweep  with 
the  data  acceptable  consideration  left  to  the  off-line  software.  After  280  spec- 
trum analyzer  sweeps  an  1120-word  buffer  is  filled  (four  words/sweepx  280 
sweeps)  and  a block  transfer  to  magnetic  tape  are  accomplished.  The  logger 
ensures  that  this  block  transfer  is  completed  before  permitting  data  acquisi- 
tion to  proceed.  Thus,  a period  of  one  spectrum  analyzer  sweep  and  two 
spectrum  analyzer  flybacks  are  occupied  with  the  block  transfer  of  data;  thus 
one  spectrum  analyzer  sweep  is  lost  every  280  sweeps.  At  the  70-Hz  rate, 
the  1120-word  buffer  is  filled  and  written  to  tape  approximately  every  four 
seconds.  The  format  of  the  7-track  magnetic  tape  for  the  peak  intensity  logger 
is  given  in  Fig.  2-16.  In  general,  each  run  (VAD  or  aircraft  fly-by),  will  gen- 
erate a distinct  number  of  1120-word  records  (one  every  four  seconds  at  70- 
Hz  rate)  with  one  end-of-file  at  the  end  of  the  run.  At  the  end  of  the  day  (or 
data  taking  period)  two  (2)  end-of-files  are  marked  at  the  end  of  the  last  run 
on  tape. 


"Data  acceptable"  condition  is  an  indication  as  to  whether  a Doppler  burst 
above  both  velocity  and  amplitude  thresholds  occurred  during  a spectrum 
analyzer  sweep. 
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Fig.  2-16  - SEL  810A  Magnetic  Tape  Format  (Peak  Intensity  Logger) 


i 


2.5  PERTINENT  LDV  PERFORMANCE  PARAMETERS 
2.5.1  Theory 

The  performance  characteristics  of  an  LDV  system  can  be  quoted  in 
terms  of  the  signal -to-noise  ratio  of  the  output  and  the  dimensions  of  the 
sensing  volume  interrogated.  Expressions  for  these  parameters  are  derivable 
from  an  electromagnetic  wave  treatment  of  the  laser  radiation  propagation  and 
interaction  with  the  aerosol  medium.  For  the  Cassegrainian  optical  system  of 
interest  here  the  field,  u,  generated  by  the  telescope  assuming  a TEM^q^  laser 
beam  mode  shape  is  given  by  (Ref.  2) 


u(t,x',  y',  z') 


Jq  <2  x)  expj-  [l  - ip  (1  - D]  x^l  xdx  , 


2 2 2 

where  p'  = x'  + y'  = radius  at  range  z ' from  the  telescope,  k a 2s-/A,  A = lasef 

wavelength,  a = exp(-2)  radius  of  the  TEM^^^  laser  mode  at  the  telescope  exit, 

N = number  of  photons  transmitted  per  second,  a^  = telescope  outer  radius, 

a.  = radius  of  telescope  inner  hole,  J ( ) = zeroth  order  Bessel  function, 

I 2 ’ o'  ' 

p = ira  /Xz,  % - T./1,  f = range  to  focus  in  the  atmosphere , (x' , z’)  are  dimensions 
in  the  plane  of  telescope,  c = speed  of  light,  and  t = time. 

For  a^  » a (laser  beam  width  considerably  smaller  than  the  telescope 
diameter,  i.e.,  the  telescope  does  not  truncate  the  laser  beam)  a simple  treat- 
ment of  the  total  LDV  system  is  possible.  Specifically  the  following  results 
apply  (see  Ref.  2 for  details); 

The  signal-to-noise  ratio  (S/N)  at  the  output  of  the  monostatic  continuous 
wave  LDV  system  is  given  by 

S/N  = nlna)  XP/4hi/B, 
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where 

r)  = the  photodetector  quantum  efficiency 

na  = the  scattering  cross  section  of  the  aerosol  at  a 
wavelength,  X (units  of  inverse  length) 

P = laser  power  transmitted 

hi/  = photon  energy  at  the  wavelength,  X,  and 

B = bandwidth  of  the  filter  being  used  to  monitor  the 
Doppler  frequency. 

Note  that  the  above  is  independent  of  the  range  to  focus  and  also  of  the 

optic  diameter.  The  physical  reason  is  that  the  sensed  volume  varies  as 

(f/R)  where  f is  the  range  to  focus  and  R the  optic  diameter.  The  number 

4 -4 

of  targets  is  therefore  of  f which  cancels  the  f dependence  of  the  return 
from  a single  target. 

Substituting  typical  parameters  in  the  S/N  equation,  namely, 
rj  = 0.3 

n<7  = 5.03  X lO*”^  ft"^  (14  mi  visibility) 

X = 10.6)1 
P = 15  W 

B = 30  kHz  (i.e.,  a filter  with  —0,5  fps  resolution) 
h|/  = photon  energy  at  10.6)1=  1-78  x 10*^^  Btu. 

we  have  an  estimate  of  45.5  dB  for  the  signal -to-noise  ratio.  This  estimate 
should  be  further  adjusted  by  no  more  than  5 dB  of  system  optical  losses. 
Such  signal-to-noise  ratios  are  routinely  observed  from  wind  returns. 

The  spatial  resolution  cell  parameters  of  interest  are  the  length  along 
the  optic  line  of  sight  and  its  diameter.  The  latter  is  always  small  for  the 
ranges  of  interest  here  (-^0.4  in)  and  will  not  be  considered  further. 
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The  along-axis  resolution  is  given  b/  2X  i /ya  for  the  untruncated 
theory  (where  a can  be  considered  to  be  the  telescope  diameter).  In  practice, 
achieving  the  condition  a is  difficult  for  a given  beam  diameter,  a,  be- 

cause optic  sizes  become  unwieldly.  Studies  have  shown  that  the  optimum 
arrangement  for  a given  optic  size  is  to  illuminate  the  telescope  such  that 
the  exp(-2)  points  of  the  intensity  pattern  coincide  with  the  extremity  of  the 
mirror.  Under  such  a condition  the  LDV  sensing  volume  length,  AL,  along 
the  optic  axis  between  the  3 dB  antenna  pattern  points  is  given  by 

Xf^ 

AL  = 4.4  — — 2 (a  = telescope  primary  radius). 

This  equation  represents  a degradation  of  a factor  of  2.2  over  the  untruncated 
case.  The  truncated  effects  on  the  signal-to-noise  ratio  {~~  3 dB  loss)  are  of 
little  concern  here  since  in  practice  under  a variety  of  weather  conditions  in 
various  parts  of  the  continental  United  States  the  signal-to-noise  ratios  have 
always  been  ample  (>  25  dB). 

For  the  one-foot  diameter  optic  size  of  interest  the  sensing  volume 
length  at  various  ranges  are; 


Range 

(ft) 

Sensing  Volume 
Length 
fftl 

100 

+0.975 

500 

+24.4 

1000 

+97.4 

2000 

+390.0 
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2.5.2  Summary  of  Pertinent  Performance  Data 
Spatial  Resolution  vs  Range  Results 

The  antenna  pattern  of  the  laser  Doppler  system  was  investigated  by 
interrogating  a target  of  known  velocity  and  position  — a rotating  disc 
at  various  ranges  from  the  LDV.  The  procedure  involved  adjusting  the 
line  of  sight  to  intersect  the  disc  at  some  radius  to  yield  a convenient  Doppler 
shift.  Range  scanning  along  this  line  of  sight  was  then  initiated.  The  disc  sur- 
face consisted  of  No.  600  grit  sandpaper  and  four  ranges  between  230  and  800  ft 
were  carefully  measured.  The  heterodyne  output  of  the  Pb-Sn-Te  detector 
was  monitored  by  the  spectrum  analyzer  adjusted  to  operate  in  the  fixed  fre- 
quency receiver  mode  and  tuned  to  the  Doppler  frequency.  As  the  secondary 
mirror  was  scanned  the  receiver  response  was  recorded  on  a strip  chart  re- 
corder along  with  the  range  command  voltage.  The  range  increment  at  which 
the  Doppler  signal  deteriorated  to  the  one-half  power  level  was  then  graphi- 
cally obtained  from  the  chart  data.  A typical  strip  chart  record  is  shown  in 
Fig.  2-17  and  the  resolution  data  calculated  from  such  profiles  are  shown  in 
Fig.  2-18  on  which  is  also  plotted  the  expected  theoretical  behavior  for  the 
nominal  system  parameters,  i.e.,  a beam  diameter  of  12  inches  at  the  telescope 
primary  (defined  at  the  exp  (-2)  intensity  points). 

S is  observed  that  the  measured  values  are  some  40%  larger  than  theo- 
retical (i.e.,  above  diffraction  limit)  which  is  indicative  of  some  system  ineffi- 
ciences.  Among  the  contributing  factors  are  respectively:  aberrations,  alignment 
difficulties,  and  beam  size  and  shape  deviations  from  nominal.  Experience 
has  shown  that  beam  shape  deformation  is  probably  the  chief  culprit  in  that 
occasionally,  when  the  laser  beam  output  is  observed  to  be  non-gaussian  (i.e., 
when  it  contains  multiple  intensity  peaks),  serious  degradation  of  the  system 
spatial  resolution  capability  occurs.  Analysis  shows  abberation  effects  to  be 
small  and  also  that  beam  size  deviations  from  nominal  are  not  overly  important 
(a  fact  not  obvious  from  the  theoretical  equation  for  the  resolution,  i.e.,  AL  = 

, since  the  degree  of  beam  truncation  by  the  telescope  primary  has  been 

2-32 


factored  out).  Greater  consideration  of  laser  beam  quality  is  merited  during 
the  design  of  future  systems. 

Sisnal-tO'Noise  Ratio  Performance 


The  signal -to -noise  ratio  performance  of  the  laser  Doppler  unit  is  moni- 
tored by  interrogating  a target  of  known  reflectance,  the  same  rotating  disc 
as  used  for  the  spatial  resolution  measurement.  The  reflectance  of  the  No.  600 
grit  sandpaper  surface  was  determined  by  comparing  LDV  system  efficiency 
while  interrogating  such  a surface  with  a similar  measurement  using  a specially 
prepared  'flowers  of  sulphur'  surface  which  is  well  documented  to  be  very  nearl> 
Liambertian  and  of  essentially  unit  reflectance  (80%)  at  the  10-micron  wavelength. 
In  this  way  the  sandpaper  surface  reflectance  was  determined  to  be  .08.  The 
surface  presented  a 45*^  angle  of  incidence  (o)  to  the  laser  beam. 

The  signal -to-noise  ratio  for  such  a scattering  configuration  is  given  by 
the  relation 


SNR  = Kj 


h I'  B L 


- detector  quantum  efficiency  = 37%  ( including  82%  detector 
window  transnilssion  loss). 

3 transmitted  laser  power  =15  watt. 

= primary  radius  (6  inches) 

= target  reflectivity  (0.08) 

= angle  of  incidence  on  target  = 45** 

= photon  energy  at  lO/i  = 1.88  x 10*^^  joules 

= spectrum  analyser  bandwidth  (30  kHs) 

3 range  to  target  (560  ft  in  typical  set  up) 

3 beam  truncation  loss  s 0.5 

3 optical  system  transmission  losses  (typically  10%). 


Substituting 


SNR  = K2  X 2.2  X 10”, 

3 83.5  dB  less  optical  transmission  losses. 


Two  way  optical  efficiency  losses  were  determined  to  be  11.2  dB  by  comparing 
power  measurement  at  the  focus  in  the  atmosphere  to  a measurement  of  the 
laser  output  power, resulting  in  a 'target'  SNR  of  72.3  dB. 

Observed,  and  what  was  considered  satisfactory  for  operational  usage 
during  this  investigation,  was  a SNR  of  60  dB  at  the  560  ft  nominal  range.  The 
additional  12  dB  or  so  additional  loss  is  due  to  such  factors  as  optical  mis- 
alignment and  mismatch  of  beams  on  the  photodectector. 

Inasmuch  as  the  quoted  optical  system  transmission  (11.2  dB)  and  align- 
ment (12.3  dB)  losses  are  excessive,  these  represent  areas  where  the  system 
could  be  improbed.  Typical  of  the  improvements  required  are  incorporation  of: 

a.  Improved  alignment  methods 

b.  Acquisition  of  improved  optical  components 

c.  Updating  the  interferometer  design,  especially  in  the 
area  of  beam  interaction  with  the  photodetector. 

With  the  system  "peaked  up"  such  that  the  signal-to-noise  ratio  obtained  while 
interrogating  the  rotating  disc  was  60  dB,  a degree  of  performance  was  ob- 
tained which  was  satisfactory  for  wind  shear  and  vortex  detection  at  both 
Huntsville  and  JFK.  Signal-to-noise  ratios  of  wind  returns  were  generally 
in  the  range  of  10-20  dB . 

2.6  DISCUSSION  OF  LIMITATIONS  OF  PRESENT  SYSTEM 

The  system  as  presently  configured  has  several  limitations  some  of 
which  were  known  before  the  investigation  was  initiated  and  some  became 
more  obvious  during  contract  performance.  They  were  not  corrected  during 
the  effort  due  to  lack  of  availability  of  operational  components,  time  and 
budget  constraints,  lack  of  the  required  technology,  etc.  The  technology  areas 
represented  by  the  limitations  should  be  given  prime  consideration  in  the  de- 
velopment of  any  future  systems.  These  areas  are  now  considered. 


> 
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Data  Rate 

The  present  system  is  capable  of  a 70-Hz  acquisition  rate  of  velocity 
samples  which  is  marginal,  for  example,  in  the  vortex  tracking  mode.  Range 
scanning  at  a 5-Hz  rate  over  a 1000  ft  range  results  in  range  'smear- 
ing' of  approximately  150  ft.  Typically  during  the  JFK  testing  the  range  scan 
limits  were  separated  by  approximately  500  ft  with  scanning  at  a rate  of  4 Hz 
resulting  in  a range  smearing  of  approximately  60  ft.  There  is  room  for  at  least 
an  order  of  magnitude  improvement  here  which  can  be  partly  achieved  by  in- 
creasing the  sweep  rate  in  the  present  processing  technique  of  using  a sweeping 
spectrum  analyzer  to  interrogate  the  Doppler  spectrum.  Increasing  this,  of 
course,  would  place  an  additional  burden  on  the  data  acquisition  system  — 
especially  the  minicomputer  currently  in  use.  The  ultimate  solution  is  a parallel 
filter  system  or  filter  bank  with  the  integration  time  chosen  to  be  sufficiently 
short  to  allow  for  the  required  data  rate,  say  1 kHz.  This  device  would  also  in- 
crease the  duty  cycle  of  the  processing  system  (Doppler  tracker)  to  essentially 
unity,  an  increase  of  two  orders  of  magnitude  over  the  present  system,  thus 
yielding  improved  signal -to -noise  ratio  performance.  The  increased  data  rate 
would  also  accommodate  an  increased  rate  of  conical  scanning.  The  ultimate 
improvements  suggested  here  would  represent  a major  addition  to  the  present 
hardware. 

Doppler  Signal  Processing 

The  duty  cycle  of  the  present  Doppler  signal  processor  (a  scanning  spec- 
trum analyzer)  is  given  by  the  ratio  of  the  spectral  region  spanned  to  the  in- 
stantaneous filter  bandwidth,  typically  2 MHz  and  30  kHz,  respectively.  This 
approximately  1,5%  duty  cycle  could  be  improved  by  the  incorporation  of 
an  integrating  filter  bank  as  alluded  to  above. 

Rate  of  VAD  Sampling 

A VAD  scan  at  a single  altitude  is  accomplished  every  5 seconds  (0.2  Hz) 
in  the  present  hardware,  which  results  in  a 40-second  time  interval  for  inter- 
rogating eight  altitudes.  For  operational  usage  this  is  too  long  since  it  can, 
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for  exa-.nple,  yield  false  estimates  of  vertical  shear  because  of  the  time  lag 
in  interrogating  the  various  altitudes.  An  increase  in  the  rate  of  scanning  is 
mechanically  possible  and  is  consistent  with  system  performance  criteria. 

A limit  to  the  conical  scan  rate  does  exist,  however,  since  the  antenna  pattern 
should  not  translate  considerably  in  the  round-trip  radiation  time.  The  limiting 
angular  rotation  rate  which  ensures  that  the  antenna  pattern  does  not  move  by 
more  than  half  its  width  during  the  round  trip  time  is  given  by 

. 1.22  Kc 
- 4 RD  0 ’ 

where  A = laser  wavelength  (10.6  microns),  c = velocity  of  light,  R = range  to 
focus,  D = optic  diameter,  9 = conical  scan  angle.  For  R = 1154  ft  (corre- 
sponding to  an  altitude  of  1000  ft),  D = 1 ft,  and  6 = 0.5  radian,  the  limit  is 
ca<  18  rad/sec  (<^  3 Hz).  A conical  scan  capability  of  up  to  3 Hz  could  there- 
for be  incorporated  in  the  system.  The  change  would  be  relatively  minor  and 
would  reduce  the  time  of  interrogating  eight  altitudes  to  less  than  10  seconds. 

System  Range  Resolution 

The  only  fundamental  limitation  of  the  present  system  (i.e.,  monostatic 
focused)  is  its  spatial  resolution  capability  at  the  longer  ranges.  Without  a 
radical  change  in  system  concept  this  can  only  be  improved  by  increasing  the 
size  of  the  transmitter  optics.  The  spatial  resolution  at  a given  range  im- 
proves as  the  square  of  the  size  of  the  optics. 

Updating  the  existing  system's  spatial  resolution  capability  would  re- 
quire incorporation  of  techniques  to  increase  the  time -bandwidth  product  of 
the  outgoing  beam  by  imposition  of  modulation  on  the  transmitted  waveform 
such  as  through  pulsed,  FM  CW,  or  phase  coding  methods.  This  ''  ould  repre- 
sent a major  modification. 

Wind  Direction  Ambiguity 

The  system  as  presently  configured  yields  a value  for  the  wind  direction 
that  contains  an  ambiguity  of  r.  This  is  because  the  system  is  homodyne  and 


I 
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not  superheterodyne  in  nature.  The  solution  to  the  problem  is  to  offset  in 
frequency  the  local  oscillator  so  that  there  is  never  any  folding  of  the  Doppler 
spectrum  about  zero  frequency.  A device  to  accomplish  this  function,  an 
acousto-optic  modulator  (Bragg  cell),  was  to  be  incorporated  in  the  system 
prior  to  the  system  evaluation  tests  but  its  suboptimal  performance  precluded 
this.  Based  on  current  research  efforts  it  appears  that  the  problems  are 
now  resolved. 

Because  of  the  mandatory  nature  of  having  a device  in  the  system  to 
resolve  the  direction  ambiguity,  the  device  is  described  here  in  considerable 
detail  along  with  some  performance  results  obtained  from  its  use.  Reference 
is  made  to  Fig.  2- 19. 

A lithium -niobate  transducer  is  bonded  to  a germanium  crystal  (cut  along 
the  axes  shown);  the  acoustic  vibrations  generated  by  the  transducer  are  ab- 
sorbed by  a lead  backing  plate  attached  to  the  opposite  side.  The  transducer 
is  excited  at  the  desired  offset  frequency.  A laser  beam  propagating  through 
the  crystal  undergoes  reflection  at  the  acoustic  wavefronts  which  results  in 
the  reflected  beam  being  upshifted  in  frequency  relative  to  the  initial  laser 
beam  frequency.  For  the  precise  geometry  shown,  i.e.,  the  angle  between  the 
acoustic  wavefronts  and  the  incident  laser  beam  in  the  crystal  is  the  Bragg 
angle,  a single  definite  beam  at  the  upshifted  frequency  is  obtained.  It  is  this 
beam  that  is  used  as  the  local  oscillator  beam  as  shown  in  the  optical  arrange- 
ment in  Fig.  2-20. 

For  design  purposes  it  is  important  to  consider  the  origin  of  the  various 
frequency  components  in  the  local  oscillator  and  receiver  optical  paths,  respec- 
tively. They  are  as  follows: 


For  the  geometry  shown. 


Acoustic  Modulator  Included 


Local  oscillator  beam 

— feedthrough  from  the  input  laser  beam  of  amplitude  a 

f^  + fg  — desired  upshifted  beam  of  amplitude  b 

f - f — undesired  downshifted  beam  of  amplitude  c due  to 

the  acoustic  reflection  from  germanium-lead  interface. 

Receiver  beam 

f^  — backreflection  from  secondary  mirror  of  amplitude  d 

f^  + f j — desired  Doppler  component  of  amplitude  e.  f . being 
the  Doppler  shift  due  to  moving  target. 

The  Spectrum  that  results  at  the  output  of  the  photodetector  is  summarized 
in  Fig.  2-21  which  also  summarizes  considerations  regarding  the  signal  level 
budget.  The  relative  magnitudes  indicated  are  generally  borne  out  in  the 
observed  spectrum  analyzer  indications  (see  Fig.  2-22a).  Satisfactory  sep- 
aration between  the  peaks  at  (f^  - f^)  (desired)  and  at  (f^  + f j)  is  achievable 
(see  Fig.  2-22b). 

The  problem  encountered  with  the  device  involved  an  insertion  loss  of 
12  do  which  precluded  the  observation  of  wind  signals.  This  was  attributed 
initially  to  transducer  bonding  deficiencies,  a problem  that  seems  to  be 
solved.  Satisfactory  wind  returns  using  the  device  have  been  observed.  The 
remaining  problem  is  a saturation  of  the  receiver  amplifier  due  to  the  strong 
feedthrough  at  f^.  A notch  filter  is  currently  being  designed  to  alleviate 
the  problem. 
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3.  COMPUTER  SOFTWARE  SYSTEM  DEVELOPMENT 

Acquisition  and  processing  of  the  LDV  signature  is  accomplished  by 
means  of  a compact  data  handling  system  developed  specifically  for  the 
Lockheed -Huntsville  LPV . The  general  elements  of  the  LDV  data  acquisi- 
tion and  data  processing  system  are  shown  in  Fig.  3-1.  The  digitized  LDV 
intensity  versus  frequency  signal  along  with  its  coordinates  in  space  are  fed 
into  the  SEL  810  minicomputer.  Reprocessing  of  the  LDV  signal  is  carried 
out  on  the  minicomputer  utilizing  on-line  computer  programs  written  in  SEL 
machine  language.  Information  from  the  SEL  810  is  stored  on  magnetic  tape 
and  is  used  as  an  input  to  the  off-line  processing  algorithms.  Off-line  process- 
ing of  the  LDV  signal  is  carried  out  on  a Univac  1108  computer  with  programs 
written  in  FORTRAN  language  and  using  card  inputs  with  information  from 
the  logs  to  supplement  the  data.  The  flow  of  data  from  the  LDV  is  sketched 
in  Fig.  3-2  showing  both  the  on-line  and  off-line  data  processing  routines. 
On-line  manipulation  of  the  data  is  carried  out  by  the  SEL  Data  Logger  pro- 
gram. The  off-line  processing  is  carried  out  by  the  VAD  and  Vortex  Track 
program.  The  final  output  consists  of  printouts,  plots,  and  vortex  track 
tapes.  A description  of  the  data  logger  and  the  VAD  and  vortex  track  pro- 
gram and  their  operational  characteristics  is  given  in  the  following  sections. 

3.1  DESCRIPTION  OF  LDV  SOFTWARE  SYSTEM 

Data  acquisition  in  the  LDV  is  carried  out  by  the  SEL  Data  Logger 
program.  The  Data  Logger  program  preprocesses  and  records  the  LDV 
signal.  A flow  chart  of  the  Data  Logger  program  is  given  in  Fig.  3-3.  For 
each  sweep  of  10,  20  or  50  millisecond  duration,  the  Data  Logger  saves  the 
maximum  amplitude  LDV  signal,  I,  and  its  corresponding  frequency, 
which  is  above  the  amplitude  and  frequency  thresholds.  The  definition  of 
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Data  Flow  from  Lockheed  LDV 


Data  Logger  Macro 


r 


I and  and  the  shape  of  the  characteristic  LDV  spectrum  is  shown  in 

Fig.  3-4.  It  can  be  seen  that  is  the  average  velocity  associated  with 

the  flow  phenomena. 

The  output  from  the  Data  Logger  program  consisting  of  as  a func- 

tion of  time  and  space  is  shown  in  Table  3-1.  Note  that  the  type  of  information 
provided  in  channels  3 and  4 depends  on  the  type  of  mode,  VAD  or  vortex  scan, 
selected.  From  the  output  of  the  Data  Logger,  the  wake  vortex  velocity  field 
or  wind  field  can  be  reconstructed  using  off-line  processing  routines. 

Final  processing  of  the  LDV  measurements  is  carried  out  by  the  VAD 
and  Vortex  Track  program  is  shown  in  Fig.  3-5.  In  this  off-line  program 
the  array  of  V^^  values  which  is  a function  of  time  and  space  is  processed 
to  yield  tlip  three-dimensional  wind  field  (VAD  mode)  or  the  aircraft  wake 
vortex  trajectories  (vortex  mode).  The  processing  of  the  VAD  measurements 
involves  the  computation  of  the  u,  v,  w wind  components  from  the  character- 
istic sinusoidal  VAD  LDV  signature  discussed  earlier  in  Section  2 and  described' 
in  more  detail  in  Appendix  B.  The  program  is  geared  to  handle  both  the  trans- 
lated and  nontranslated  LDV  signal.  However,  during  the  course  of  this  research 
effort  all  of  the  data  acquisition  and  data  processing  was  done  in  the  non -translate 
mode  using  an  algorithm  which  processes  the  peak  signals  from  the  VAD  spec- 
trum. The  final  output  Is  a plot  (and  printout)  of  the  u,v,  and  w velocity  com- 
ponents as  a function  of  altitude  and  time.  The  peak  algorithm  was  developed 
under  the  contract  to  process  the  JFK  VAD  measurements.  Subsequently,  the 
capability  was  added  to  the  processing  software  to  consider  spectral  process- 
ing for  the  winds  using  the  inverted  spectrum  and  also  using  a sine  curve  fit. 

For  the  sake  of  completeness,  the  details  of  the  full  VAD  processing 
are  given  in  Appendix  B and  are  outlined  in  this  report. 
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Fig.  3-4  - Typical  Vortex  Spectrum 
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OUTPUT  FROM  DATA  LOGGER  PROGRAM 


NOT  USED;  ALWAYS  OFF 


NORMAL  AND  CONICAL  SCAN 


WORD  FORMAT:  NORMAL  SCAN 


max 


Table  3-1  - (Continued) 


CHANNEL  4 (OCTAL  63) 


-64  (BIT  1 OF  CHANNEL  1 IDENTIFIES  WHICH) 


^4 

I 

I 

1 

I 

! 

When  the  LDV  system  is  operating  in  the  finger  scan  vortex  track  mode, 
the  VAD  and  Vortex  Track  program  flow-charted  in  Fig.  3-5  calculates  the  alti- 
tude and  lateral  location  of  the  port  and  starboard  aircraft  vortices  as  a function 

of  time  from  the  array  of  V values.  The  vortex  location  algorithm  cal- 
' ms  ” 

culates  and  sorts  the  line -of-sight  velocity  as  a function  of  elevation  angle 

and  time,  selects  the  highest  line-of-sight  velocities  or  "hits"  observed  during 

each  frame,  computes  the  vortex  location  from  the  centroid  of  the  high  velocity 

"hits,"  and  plots  the  track  of  the  vortex  centroid.  In  addition,  scatter  plots 

and  printouts  of  the  individual  hits  are  generated.  A detailed  description  of 

the  vortex  tracker  subroutines  is  given  in  Appendix  B and  the  basic  sorting 

criteria  involving  V is  described  in  Ref. 4. 
ms 

3.2  OPERATION  OF  LDV  SOFTWARE  SYSTEM 

Operation  of  the  Lockheed -Huntsville  LDV  involves  initialization  of  the 
SEL  Data  Logger  program  and  the  recording  of  the  LDV  signature.  After  the 
measurements  have  been  recorded  at  the  proper  threshold  settings,  the  VAD 
and  Vortex  Track  program  is  used  to  process  the  data  into  VAD  or  vortex 
track  plots.  It  is  useful  to  examine  the  data  processing  operations  involved 
in  determining  the  wake  vortex  trajectories  and  the  three-dimensional  wind 
velocity  field  from  the  LDV  measurements  In  terms  of  on-line  preprocessing 
and  off-line  processing. 

3.2.1  On-Line  Data  Processing 

The  conical  scan  VAD  and  finger  scan  vortex  track  measurements  are 
preprocessed  and  recorded  by  the  SEL  computer.  The  input  from  the  SEL 
computer  consists  of  the  basic  signal  as  well  as  additional  test  param- 

eters which  are  listed  in  Table  3-2.  The  data  recorded  by  the  SEL  is  in  a 
different  format  depending  on  which  mode  (VAD  or  vortex  track)  the  system 
is  operating  in.  The  SEL  records  the  basic  LDV  signal  with  the  Data  Logger 
program. 
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Table  3-2 

INPUT  FROM  SEL  COMPUTER 


Note;  Each  vortex  and  VAD  run  are  in  a separate  file.  The  following  data 
are  provided  for  each  spectrum  analyzer  sweep. 


1, 

Spectral  Sample  Count  Across  Spectrum  Analyzer  Sweep, 
Corresponding  to  V 

2. 

Amplitude  at  the  Above  Point 

3. 

Data  Acceptable  Flag 

4- 

Flag  for  Spectrum  Analyzer  Sweep  Speed 

5. 

Flag  for  Translator 

6. 

Flag  for  Positive  or  Negative  Frequency  (used  only  when 
translator  is  used) 

7. 

Flag  for  V or 

* peak  ms 

8. 

Flag  for  Conical  Scan  or  Normal  Scan 

For  Conical  Scan 

9. 

Height  above  Van 

10. 

Flag  for  Azimuth  Switch 

11. 

Cone  Angle 

! For  Finger  Scan 

9. 

Range 

10. 

Elevation  Angle 

A dump  of  a sample  output  tape  from  the  Data  Logger  program  operating 
in  the  vortex  track  mode  is  shown  in  Fig.  3-6,  each  row  corresponding  to  infor- 
mation recorded  for  each  spectrum  analyzer  sweep.  The  information  is  sep- 
arated into  12  columns  in  the  printout  with  interpretation  of  the  various  columns 
as  follows; 


Column 

1 (PEAK/MAX) 


2 (9  = 0) 

3 (NOR/VAD) 

4 (SWEEP  SPEED) 

5 (DIG  TRCK) 

6 (DATA  ACCEPT) 

7 (+/-) 

8 (TRANS/NON - 
TRANS) 

9 (SPECTRUM 
INTENSITY) 

10  (RANGE) 

11  (LTRNC  TRCK) 

12  (ELEVATION 


Interpretation 

An  indication  as  to  whether  the  peak  Doppler  fre- 
quency, fpeak>  maximum  Doppler  frequency 
above  amplitude  or  frequency  thresholds  (f  ) 
algorithm  is  utilized  (see  Fig.  3-3).  (The 
algorithm  was  used  at  all  times  during  this  investi- 
gation.) 

Inoperative . 

An  indication  as  to  whether  the  system  is  operating 
in  the  normal  (vortex  track)  (1)  or  VAD  (0)  mode. 


Sweep  speed  of  the  spectrum  analyzer  trace  (1 
indicates  10  msec/cm). 


The  computer  calculated  estimate  of 
of  full  scale). 


peak 


(percent 


An  indication  that  there  was/was  not  an  output  above 
frequency  amplitude  thresholds  during  a sweep. 


An  indication  as  to  the  sense  of  the  Doppler  shift 
(i.e.,  target  moving  toward  or  away  from  trans- 
ceiver). Toward  = +,  Away  = - . 


An  indication  as  to  whether  or  not  a frequency 
translator  was  incorporated.  (During  this  investi- 
gation it  was  not  incorporated.)  No  = 0,  Yes  = 1. 

Peak  amplitude  of  the  Doppler  spectrum  in  region 
above  a frequency  threshold  (arbitrary  units). 

Range  to  focus  of  laser  Doppler  system  (ft) 

On-line  frequency  tracker  estimate  of  fpeak  (should 
be  identical  to  column  5 with  time  lag  oT  one  sample 
except  when  data  are  not  acceptable). 

Instantaneous  scan  elevation  angle  in  degrees. 
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Fig.  3- 6 - Dump  of  Sample  Output  Tape  from  Data  Logger  with  System  Operating  in  Vortex  Track  Mode 


A dump  of  a sample  output  tape  from  the  Data  Logger  operating  in  the 
VAD  mode  is  shown  in  Fig.  3-7,  each  row  of  data  corresponding  to  informa- 
tion recorded  for  each  spectrum  anal/zer  sweep.  The  information,  packed 
into  four  computer  words  on  tape,  is  separated  into  13  columns  in  the  printout 
with  interpretation  of  the  various  columns  being  as  follows; 


Column 

1 (PEAK/MAX) 


2 (9  = 0) 

3 (NOR /VAD) 

4 (SWEEP  SPEED) 

5 (DIG  TRCK) 

6 (DATA  ACCEPT) 

7 (+/-) 

8 (TRANS/NON- 
TRANS) 

9 (SPECTRUM 
INTENSITY) 

10  (NBR  ROTN) 

11  (ALT  (M)) 

12  (LTRNC  TRCK) 

13  (CONE  ANGLE) 


Interpretation 

An  indication  as  to  whether  the  peak  Doppler  fre- 
quency (fpeak^  mciximum  Doppler  frequency 
above  amplitude  and  frequency  thresholds  (f  ) 
algorithm  is  utilized  (see  Fig. 3-3).  (The  f^g  ^ 
algorithm  was  used  at  all  times  during  this 
investigation.) 

A conical  scan  azimuth  reference  which  is  nonzero 
when  the  reference  switch  is  activated. 


An  indication  as  to  whether  the  system  is  operating 
in  the  normal  (vortex  track)  (1)  or  VAD  (0)  modes. 


Sweep  speed  of  the  spectrum  analyzer  trace  in 
msec/cm. 


The  computer  calculated  estimate  of  f 


of  full  scale). 


peak 


(percent 


An  indication  that  there  was/was  not  an  output 
above  frequency  and  amplitude  during  a sweep. 

An  indication  as  to  the  sense  of  the  Doppler  shift 
(i.e.,  target  moving  toward  or  away  from  trans- 
ceiver). Toward  = +,  Away  = - . 


An  indication  as  to  whether  or  not  a frequency 
translator  was  incorporated.  During  this  investi- 
gation it  was  not  incorporated.  No  = 0,  Yes  = 1. 

Peak  amplitude  of  the  Doppler  spectrum  in  region 
above  a frequency  threshold. 

Number  of  successive  VAD  scans  for  a particular 
altitude . 


Altitude  of  VAD  for  particular  sweep 

On-line  frequency  tracker  estimate  of  fpeak  (should 
be  identical  to  column  5 with  time  lag  oi  one  sample 
except  when  data  is  not  acceptable). 

Half-angle  of  VAD  cone  in  degrees. 
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Fig.  3-  7 - Dump  of  Sample  Output  Tape  from  Data  Logger  with  System  Operating  in  VAD  Mode 


3.2.2  Off-Line  Data  Processing 

Based  on  the  array  of  line-of-sight  velocities  recorded  by  the  SEL,  the  » 

VAD  and  Vortex  Track  program  computes  the  wind  field  and/or  wake  vortex 
trajectories.  A sunnmary  of  the  basic  off-line  data  processing  technique  is 
given  in  Table  3-3.  The  data  processing  steps  described  in  Table  3-2  form 
the  fundamental  framework  for  the  VAD  and  Vortex  Track  program.  The  VAD 
and  Vortex  Track  program  is  a flexible  and  comprehensive  software  package 
which  is  the  off-line  data  processing  element  of  the  LDV  software  system. 


Table  3-3 

SUMMARY  OF  OFF-LINE  DATA  PROCESSING  TECHNIQUE 


DESCRIPTION  OF  VAD  PROCESSING  TECHNIQUE 

1. 

.Save  line-of-sight  velocities  for  one  rotation  of  scanner. 

2. 

If  two  or  more  rotations  at  same  altitude,  average  with  first 
rotation. 

3. 

Assign  azimuth  angle  to  each  point  (assuming  constant  rotation 
rate).  Plot  velocity  versus  angle. 

4. 

Edit  points . 

5. 

Apply  moving  average  if  wanted. 

6. 

Plot  velocity  versus  angle. 

PEAK  ALGORITHM  TECHNIQUE 

7. 

Pick  two  peak  velocity  points,  V_  and  Vp  , that  occur  a 

^ 1 

minimum  of  90  deg  apart. 

8. 

Compute  horizontal  velocity,  Vj^ 

Vp  + Vp 

V - 12 

h , . (cone  angle) 

^ Sin  i ^ 

9. 

Compute  horizontal  angle  with  help  of  estimated  wind  direction. 

10. 

Compute  vertical  wind  velocity,  V^ 

V - V 

p p 

I 2 

* 2 cos  l5°g«^ngiej. 

11. 

Derectify  VAD  signal  if  no  translator  is  present. 
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Table  3-3  (Continued) 


FOURIER  COEFFICIENT  TECHNIQUE 

12. 

Compute  Fourier  coefficients  (program  computes  first  four  coef- 
ficients). 

= -f  E (^n  + ‘’n 

n=:l 

= i 'tjr.  = 1 E 

where  L = number  of 
points  in  VAD  sweep. 

13. 

Compute  correction  to  fundamental  harmonic  due  to  frequency- 
cutoff  (Ref.  6). 

Correction  = ^ — ; k where  Z = No.  Zeros/No.  Points 

1 . z + ^ 

TT 

Correction  Value  = (Correction)  * (Calculated  Fundamental) 

14. 

Compute  vertical  wind  correction.  (Ref.  6) 

Correction  = ^ 

15. 

Compute  horizontal  velocity. 

• 

„ _ 

h . .cone  angle. 

8tn( 2 ® — ) 

I 
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Table  3-3  (Continued) 


16. 


17. 


Compute  horizontal  angle. 

AnglCj^  = Atan(bj/aj) 

Compute  vertical  wind  velocity 


-a 


V = 


V ■ 2 co8(^?"Y^^^) 


SINE 


CURVE  FIT  TECHNIQUE 


19. 


Find  least  squares  curve  fit  for  a sine  wave  to  the  data 

L 


20. 


Minimum  = (Y.  - 

1=1 

where 

is  line-of-sight  velocity  at  point 
9^  is  azimuth  at  point 

C,  A and  B are  coefficients  to  be  solved  for 
Compute  horizontal  velocity. 


C - A cosSj  - B ain9jl 


Va^  + B^ 


. .cone  ansle. 
sin( 2 — 


21. 


22. 


Compute  horizontal  angle. 

Angle^  = Atan(B/A) 

Compute  vertical  wind  velocity 

V a z£ 

V __  .cone  angle, 
cos( 2 — ® — ) ' 
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AFTER  THE  VELOCITIES  HAVE  BEEN  COMPUTED  FOR  ALL  ALTITUDES: 

1.  Compute  a least  squares  fit  for  the  horizontal  speed  to  a power 
law  curve 

V = (H/H^)P 


Compute  a least  squares  fit  for  the  horizontal  directions  to  a 
polynomial. 

Plot  horizontal  speed  versus  height. 

Plot  horizontal  direction  versus  height. 

5.  Plot  crosswind  velocity  versus  height. 


DESCRIPTION  OF  VORTEX  TRACKING  ALGORITHM 


1.  Frequency,  intensity,  range,  elevation  angle,  and  time  are  saved 
for  each  point  in  a frame. 

2.  The  X,  y position  is  computed  with  adjustment  for  the  wind. 

3.  The  line-of-sight  airspeed  is  computed.  If  the  translator  is  used, 
the  crosswind  speed  is  removed. 

4.  The  points  are  sorted  with  respect  to  velocity  and  also  sorted  with 
respect  to  intensity. 

5.  Depending  on  algorithm  chosen,  the  highest  velocity  point  or  the 
highest  intensity  point  is  picked. 

6.  A correlation  circle  is  drawn  around  the  above  point.  Only  points 
within  this  circle  will  be  considered  for  determining  the  vortex. 

7.  The  number  of  points  within  the  correlation  circle  are  counted 
and  compared  with  NPSUF.  If  Insufficient,  throw  out  the 
chosen  point  and  return  to  item  5. 

8.  The  fraction  of  points,  with  a velocity  or  intensity  greater  than 
APERCT  of  the  maximum  point,  is  computed.  If  this  fraction  is 
less  than  BPERCT,  throw  out  the  chosen  point  and  return  to  item  5. 

9.  The  velocity  at  each  point  is  multiplied  by  its  intensity,  to  be  used 
as  a weighting  factor  in  determining  the  centroid  of  the  points. 

This  centroid  is  declared  to  be  the  vortex  loc..tion.  The  relative 
weighting  of  intensity  versus  velocity  for  computing  the  centroid 
may  be  modified  by  the  use  of  NOtSEF  (the  noise  floor)  and  ADJI 
(intensity  adjustment  or  fraction  of  noise  floor  added  to  total 
intensity). 
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Table  3-3  (Concluded) 


10.  An  excluding  circle,  with  a radius  of  EPERCT  multiplied  by  the 
correlation  radius,  is  drawn  around  the  above  vortex. 

11.  The  highest  velocity  point  or  the  highest  intensity  point  outside 
the  excluding  circle  is  chosen.  (In  the  case  of  the  velocity  algo- 
rithm, no  point  with  less  than  half  the  velocity  of  the  highest 
velocity  of  vortex  1 will  be  considered.) 

12.  A correlation  circle  is  drawn  around  the  above  point  and  only  points 
within  this  circle  and  outside  of  a correlation  circle  around  vortex  1 
will  be  considered  for  determining  vortex  2. 

13.  The  number  of  points  within  the  area  described  above  is  counted  and 
compared  with  NPSUF  and  CPERCT  multiplied  by  (number  of  points 

I used  in  determining  vortex  1).  If  insufficient,  throw  out  the  chosen 

point  and  return  to  item  11. 

14.  The  fraction  of  points  with  a velocity  or  intensity  greater  than 
APERCT  of  the  maximum  point  is  computed.  If  this  fraction  is 

1 less  than  BPERCT  throw  out  the  chosen  point  and  return  to  item  11. 

15.  The  velocity  at  each  point  is  multiplied  by  its  intensity  to  be  used 
as  a weighting  factor  in  determining  the  centroid  of  the  points. 

This  is  declared  to  be  the  vortex  2 location. 

I 16.  The  vortex  with  a greater  Y position  is  declared  to  be  the  star- 
board vortex  with  the  other  vortex  being  port.  If  only  one  vortex 
is  found  it  is  declared  unknown. 

17.  A scatter  plot  showing  the  magnitude  and  location  of  each  velocity 
point  is  made  for  each  frame  on  the  printer. 

18.  SC  4020  plots  and  printer  plots  are  made  of  the  following 

a.  Vortex  height  versus  horizontal  position 

b.  Vortex  height  versus  time 

c.  Vortex  horizontal  position  versus  time. 


From  the  type  of  information  shown  earlier  in  Figs.  3-6  and  3-7,  the 
VAD  and  Vortex  Track  program  is  used  to  reconstruct  the  ambient  wind  field 
or  the  aircraft  wake  vortex  trajectories.  To  describe  the  operation  of  the 
software,  it  is  useful  to  present  sample  runs  and  to  analyze  the  results.  The 
list  of  inp<it  parameters  to  the  off-line  analysis  program  and  their  definitions 
are  provided  in  Fig.  3-8.  A sample  VAD  run  from  the  VAD  and  Vortex  Track 
program  showing  the  resulting  printout  and  plots  is  given  in  Figs,  3-9  and 
3-10,  respectively. 
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TiTLE  CaRU  3z  COLUMNS  Op  cOMmEuT 

offline  OATa  reduction  program  pOK  LMSc  LASEk  Van 
handles  vao  aino  Scan  and  normal  vortex  arc  scan 

program  input 


namelist  (oATaI 


IGROUP  NUMBER  OP  GROUPS  qF  OAT*  TO  BE  SxiPPED 

isfilEi2,io)  skip  files  isfiLe(1,k)  Through  isfileiz.ki 

NRUN  NUMBER  Up  PLYBYS  aND  VaD  RUNS  TO  pROCESS 

(Includes  skipped  piles  but  not  Skipped  groups) 

(IF  NKUN«0  program  IS  TERMINATED) 

zlaser  height  above  ground  of  Laser  m mirRor  for  vortex  sca 

ZLASCN  height  above  ground  of  Conical  scan  mirror 

intvel  flag 

INTVEL  - ( velocity  ORIenTEO  VorTex  DETERMINATION 

inTvel  ■ 2 intensity  Oriented  voRtex  determination 
NPsuF  sufficient  Number  of  points  to  determine  vortex  position 

aperct  fraction  of  The  maximum  peax  y 

BPERCT  fraction  of  points  of  which  the  0 IS  at  least 

APERCT  fraction  of  The  Maximum  u 
(8  IS  velocity  or  intensity  as  determined  by  INTvEl 
cperct  fraction  of  points  Used  in  vortex  one  required  for 

VORTEX  two 

RPERCT  fraction  OF  AIRCraFT  WINg  SPAN  USED  FOR 

correlation  HAOluS 

represented  aircraft  aR£  8707,  8727,  b737,  b7R7, 
oca,  DC9,  ocic>  l>dii,  C5A,  cvaeo 
epekct  fraction  of  correlation  radius  from  vortex  one  for 

excluding  Initial  point  of  vortex  two 

noisef  noise  floor 

AuJI  intensity  aOjusTmeNI 

(FRACTION  OF  noise  FlOOR  ADDED  TO  TOTAL  INTENSITY) 

ANGSW  A^MUTH  of  MIrROR  wHEn  switch  is  activated 

wangle  estimate  of  wind  angle 

(direction  wind  is  f Rom ) 

(USEO  when  no  translator  is  USED) 
windhp  wind  height  to  Be  Plotted 

(IF  zero  wind  will  “e  plotted  to  height  of  conical 

SCAN! 

LFLIP  FLAG  (aPPLTS  ONLY  FOR  NOnTRANSLATE  MODE) 

lflip  • c USES  peaks.  Does  not  flIp  data 


Fig.  3-8  - Input  Parameters  to  VAD  and  Vortex  Track  Program 
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Fig.  3-9  - Sample  of  VAD  and  Vortex  Track  Program  Operating  in  VAD  Mode 
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Fig. 3-9  - (Concluded) 


Fig.  3-10  - Sample  of  VAD  and  Vortex  Track  Program  Operating 
in  VAO  Mode 

a.  Basic  Signal 
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Fig.  3-10  (Concluded) 

c.  Computed  Flip  of  Basic  Signal 


The  plots  corresponding  to  the  sample  printout  shown  in  Fig.  3-9  are 
presented  in  Fig.  3-10.  The  plots  include  the  observed  line-of-sight  velocity 
as  a function  of  azimuth  angle,  the  signature  filtered  by  a 21 -point  moving  aver- 
age, and  the  signal  derectlfied.  For  the  latter  two  cases,  the  signal  is  edited. 

The  edit  criterion  is  applied  to  all  of  the  points  whose  magnitude  is  above  the 
specified  velocity  threshold.  The  edit  program  compares  successive  points 
in  the  velocity  versus  horizontal  scan  angle  distribution  and  deletes  those 
points  where  the  magnitude  of  the  velocity  differs  by  more  than  a specified 
fraction  (20%  for  the  JFK  tests)  from  the  adjacent  values.  The  purposes  of 
this  edit  criteria  is  to  filter  out  high  frequency  turbulence  and  to  isolate  the 
fundamental  mean  flow  components. 

The  printout  from  the  off-line  program,  illustrated  in  Fig.  3-9,  con- 
sists of  the  list  of  input  parameters,  test  parameters,  and  VAD  measurements. 
The  VAD  measurements  show  the  wind  speed  and  wind  direction  as  a function 
of  altitude  and  time.  In  the  sample  output  the  wind  speed  and  direction  at 
each  altitude  (and  time  interval)  is  computed  three  different  ways  according 
to:  (1)  the  peak  algorithm,  (2)  the  spectral  analysis,  and  (3)  the  sine  curve  fit. 

The  winds  are  computed  with  the  peak  algorithm  for  the  two  highest 
line-of-sight  velocities  observed  in  the  VAD  (resulting  in  wind  speed  and 
direction  1 and  2,  respectively,  in  Fig.  3-9).  The  average  of  wind  speed 
and  direction  1 and  2 is  also  given  by  the  peak  algorithm  and  is  assumed  to 
be  more  representative  of  the  actual  phenomenon.  In  the  printout  the  re- 
versed wind  direction  is  also  given  since  a 180  degree  ambiguity  is  present 
in  the  non -translated  LDV  mode  which  is  resolved  by  specifying  the  wind 
quadrant  as  an  input.  The  lower  and  upper  angle  refer  to  the  angle  over 
which  the  VAD  signal  is  flipped  or  derectlfied  (see  Fig.  3-  10). 

The  spectral  analysis  program  computes  the  wind  speed  and  direction 
based  on  the  first  harmonic  of  the  VAD  signal  as  shown  in  the  results  in  Fig. 
2-10.  In  addition,  the  seconu  and  third  harmonic  contributions  are  computed 
and  printed  out.  The  number  of  points  observed  in  the  VAD  scan  (the  number 
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of  hits)  and  the  number  of  zeros  are  tabulated  and  printed  out.  The  zero  factor 
is  a scale  correction  factor  to  adjust  the  wind  speed  of  the  fundamental  biased 
by  the  frequency  threshold.  As  shown  in  Fig.  3- 10,  when  the  VAD  sigral  is 
derectified  gaps  occur  at  the  crossover  regions  due  to  the  frequency  (velocity) 
threshold.  The  zero  factor  is  computed  to  adjust  the  wind  speed  based  on  the 
number  of  zeros  in  the  VAD  scan  (i.e.,  in  Fig.  3-9  the  wind  speed  computed 
from  the  first  harmonic  was  multiplied  by  1.5148  to  derive  the  actual  wind 
speed).  The  RAW  VERTICAL  DATA  refers  to  the  measured  vertical  velocity 
without  the  frequency  cutoff  correction,  and  the  VERTICAL  WIND  SPEED  in 
refers  to  the  corrected  value  as  shown  in  Fig.  3-9. 

The  sine  curve  fit  algorithm  computes  the  wind  speed  and  di:.ection 
making  a least  square  fit  to  the  VAD  signature  with  a sine  wave  and  the  re- 
sults are  shown  in  Fig.  3-9.  The  standard  deviation  of  the  Une-of-sight 
velocity  about  the  sine  wave  is  also  computed  and  given. 

Data  processing  of  the  finger  scan  LDV  measurements  is  carried  out 
by  the  vortex  track  routines  in  the  VAD  and  Vortex  Track  program.  The 
operation  of  the  vortex  tracking  software  was  summarized  earlier  in  Table 
3-1,  and  is  illustrated  by  the  sample  output  shown  in  Fig.  3-  11.  Typical 
printout  from  the  program  includes  the  number  of  saved  hits  (IP)  and  the 
observed  signal  frequency  (IFREQ  = V^^)  which  is  tabulated  as  a function 
of  time,  range,  and  elevation  angle  for  each  scan.  The  line-of-sight  velocity 
(SPEED  in  ft/sec)  and  lateral  and  vertical  coordinates  (YP  and  ZP  in  feet) 
are  computed  from  the  above  parameters.  The  points  are  sorted  according 
to  their  intensity,  selected  by  the  tracking  algorithm  for  the  determination 
of  the  centroid  of  the  port,  starboard,  or  undefined  vortex  (YC,  ZC).  For 
each  elevation  sweep  a scatter  plot  of  the  intensity  of  all  of  the  hits  is  given 
labeled  A through  O and  a list  of  points  used  in  determining  the  centroids  is 
given.  In  the  scatter  plots  the  Z's  represent  the  points  about  which  the  corre- 
lation circle  are  drawn  and  the  P's  and  S's  are  the  centroids  of  the  correlation 
circles  and  denote  the  location  of  the  port  and  starboard  vortices,  respectively. 
The  line-of-sight  velocity  is  processed  and  displayed  in  this  manner  for  each 
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successive  sweep  until  the  end  of  the  data  is  reached.  At  this  point,  the  list 
of  vortex  centroids  as  a function  of  time  is  printed  out  and  plotted.  The  wake  ♦ 

vortex  trajectories  are  also  plotted  on  a SC  4020  plotter.  The  wake  vortex 
plots  generated  from  the  JFK  tests  are  discussed  in  Section  5 and  the  measure- 
ments are  presented  in  Appendix  A. 

The  manner  in  which  the  wake  vortex  measurements  were  processed 
can  be  summarized  as  follows.  The  frequencies  and  amplitudes  associated 
with  the  laser  Doppler  signal  were  sampled  at  fixed  intervals.  The  spectrum 
was  recorded  if  it  was  above  the  frequency  and  amplitude  threshold  settings 
(Fig.  3-4).  The  amplitude  and  frequency  threshold  settings  for  the  JFK  tests 
are  given  in  the  log  sheets  in  Appendix  C.  From  the  array  of  recorded  fre- 
quency and  intensity  points,  the  line-of-sight  velocity  field  was  computed  and 
the  vortex  parameters  including  location  and  velocity  distribution  were  deter- 
mined. 

To  compute  the  wake  vortex  transport  and  decay  characteristics  from 
the  low-speed  data  line-of-sight  velocity  distributions,  the  JFK  measurements 
were  analyzed  using  the  VAD  and  Vortex  Track  Program  described  earlier. 

Based  on  previous  experience  with  the  program,  the  follovring  parameters 
were  selected  for  the  analysis  of  the  JFK  data: 


ii 


INTVEL  = 2.0 


NPSUF  = 4.0 
APERCT  = 0.1 
BPERCT  = 0.1 


CPERCT  = 0.5 

RPERCT  = 0.314 


FUg 

INTVEL  = 1 Velocity  oriented  vortex  determination 
INTVEL  = 2 Intensity  oriented  vortex  determination 

Sufficient  number  of  points  to  determine  vortex 
position 

Fraction  of  the  maximum  peak  velocity  or  intensity 
points 

Fraction  of  points  within  the  correlation  circle 
where  Q is  at  least  APERCT  fraction  of  the 
maximum  Q (Q  is  velocity  or  intensity  as  deter- 
mined by  INTVEL) 

Fraction  of  number  of  points  in  correlation  circle 
used  for  determining  vortex  1 required  for  deter- 
mination of  vortex  2 

Fraction  of  aircraft  wing  span  used  for  correlation 
radius 
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EPERCT  = 1.5 

NOISEF  = 0.0 
ADJI  = 0.0 


Fraction  of  correlation  radius  from  vortex  1 for 
excluding  initial  point  of  vortex  2 

Noise  floor 

Intensity  adjustment  (Fraction  of  noise  floor 
added  to  total  intensity) 
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4.  FIELD  TESTS  OF  LDV  SYSTEM 

To  evaluate  the  capability  of  the  Lockheed -Huntsville  LDV  to  measure 
the  wind,  wind  shear,  and  wake  vortex  phenomena  in  terminal  areas;  the  sys- 
tem was  deployed  at  the  Huntsville  Jetplex,  Huntsville,  Alabama,  and  at 
Kennedy  International  Airport,  Jamaica,  New  York. 

The  LDV  was  deployed  and  tested  at  the  Huntsville  Jetplex  in  October 
1975  to  check  out  the  basic  operation  of  the  system.  The  location  of  the  LDV  at 
the  Jetplex  was  as  shown  in  Fig.  4-1.  On  15  October  1975  wake  vortex  measure- 
ments were  made  behind  an  L-1011  aircraft  and  the  following  day,  measure- 
ments were  made  behind  a FAA  CV-880  aircraft  conducting  touch  and  go 
landings  and  approaches.  The  wake  measurements  at  the  Huntsville  Jetplex 
were  intended  primarily  to  test  the  operation  of  the  system.  The  ability  of 
the  LDV  system  to  identify  and  track  wake  vortices  was  demonstrated  at  the 
Jetplex  tests.  A sample  wake  vortex  trajectory  measured  by  the  LDV  at  the 
Jetplex  is  shown  in  Fig.  4-2.  The  vortex  track  in  Fig. 4-2  shows  the  lateral 
motion  of  the  L-1011  wake  vortices.  Similar  wake  measurements  were  ob- 
tained at  Huntsville  for  the  CV-880  aircraft.  The  ability  of  ttte  LDV  system 
to  monitor  wake  vortices  at  an  airport  for  extended  periods  of  time  under 
various  weather  conditions  was  demonstrated  during  the  JFK  tests. 

Wake  vortex  and  wind  profile  nMasurements  were  taken  with  the 
Lockheed-Huntsville  LDV  at  Kennedy  International  Airport  (JFK)  in  two 
test  sequences  from  30  October  1975  to  19  November  1975  and  from  13 
January  to  31  January  1976.  (Between  those  tests,  the  LDV  was  transported 
to  Rosamond.  California,  to  participate  in  the  DOT /NASA  B-747  wake  vortex 
decay  studios.)  In  the  14  days  of  operation  at  JFK  in  November  1975.  13 
tapes  of  vortex.  VAD,  and  range  scan  data  were  generated.  During  the 
second  test  series  in  January  1975  anotiier  13  tapes  of  VAD  and  wake  vortex 
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Fig,  4-1  - Field  Site  Location  for  Mimtsville  Jetplex  LDV  Teats 
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measurements  were  obtained.  A total  of  479  landing  operations  on  Runway 
31R  involving  B-747,  B-727,  B-707,  DC-10,  DC-9,  DC-8  and  L-1011  aircraft 
were  monitored  by  the  LDV  system  at  JFK.  The  results  of  the  JFK  tests 
Are  discussed  in  Section  5 and  the  measurements  and  data  logs  are  given  in 
Appendixes  A and  C,  respectively.  The  location  and  orienUtion  of  the  LDV 
with  respect  to  the  runway  at  JFK  is  sketched  in  Fig.  4-3. 


5.  RESULTS  OF  WAKE  VORTEX  MEASUREMENTS 

To  demonstrate  the  operational  capability  of  the  LDV  for  monitoring 
wake  vortices  at  airports,  wake  vortex  measurements  were  carried  out  with 
the  LDV  system  for  479  landing  operations  on  runway  31R  at  Kennedy  Inter- 
national Airport  (JFK).  From  the  logged  data,  460  flybys  have  been  processed 
and  the  presence  of  the  wake  vortex  was  clearly  detected  with  the  system  for 
379  flybys  (i.e.,  one  or  more  vortex  centroids  were  computed)  which  includes 
82%  of  the  processed  data.  For  the  remaining  18%  of  the  data  (81  cases)  no 
clear  identification  of  the  wake  vortex  pair  was  made.  The  lack  of  vortex 
identification  is  attributed  to;  (1)  light  aircraft  types  (such  as  DC-9  and  others) 
where  the  vortex  strength  is  too  small  for  the  LDV  to  detect  them;  (2)  operator 
error  in  initializing  the  start  of  the  scan  or  in  setting  the  frequency  and  ampli- 
tude thresholds;  and  (3)  windy  and/or  gusty  conditions  where  the  wake  vortex 
is  rapidly  dispersed  or  transported  out  of  the  field  of  view.  A summary  of  the 
wake  vortex  measurements  is  listed  in  Table  5-1  including  the  date,  the  number 
of  flybys  logged,  the  number  of  flybys  processed,  the  number  of  observed  wake 
trajectories,  and  pertinent  comments.  The  individual  wake  vortex  trajectories 
are  given  in  Appendix  A for  those  flybys  out  of  the  total  of  379  where  the 
vortex  tracks  were  within  the  field  of  view  of  the  LDV  system  (i.e.,  six  or 
more  vortex  centroids  were  located).  The  general  trends  are  described  below. 

Sample  wake  vortex  trajectories  recorded  with  the  LDV  at  JFK  are  pre- 
sented in  Figs.  5-1  through  5-10.  These  flybys  have  been  selected  since  they 
Illustrate  the  wake  vortex  measuring  capability  of  the  systems  in  terms  of: 

a.  Discriminating  the  port  and  starboard  vortices  from  the  overall 
flowfield, 

b.  Tracking  the  laternl  and  vertical  location  of  the  vortex  wake, 
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Table  5-1 

SUMMARY  OF  WAKE  VORTEX  MEASUREMENTS  CONDUCTED 
AT  KENNEDY  INTERNATIONAL  AIRPORT  WITH 
LOCKHEED -HUNTSVILLE  LDV 
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rig.  S*1  • Wak«  VortM  Trajectory  of  B-747  Aircraft  — Date, 
10/31/7S;  Time.  16t51t30 
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Fig.  5-2  - Wake  Vortex  Trajectory  of  B-747  Aircraft  - Oate, 
10/31/75;  Time,  17:15:00 
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Fig,  5-5  - Wake  Vortex  Trajectory  of  B-747  Aircraft  — Date  11/18/75 
Time.  17:2:00 
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6 - Wake  Vortex  Trajectory  of  B-747  Aircraft  — Date,  1/14/76; 
Time,  16:47:00 
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Fig.  5-8  - Wake  Vortex  Trajectory  of  B-747  Aircraft  — Date,  1/30/76; 
Time,  15:15:00 
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Fig.  5-9  - Wake  Vortex  Trajectory  of  B-747  Aircraft  - Date,  l/30/7<> 
Time,  16:10:00 
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Fig, 5-10  - Wake  Vortex  Trajectory  of  B-747  Aircraft  — Date  1/30/76; 
Time,  17:00:00 
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c.  Indicating  transport  of  the  vortex  pair,  and 

d.  Indicating  decay  of  the  vortex  peak  velocity. 

The  port  and  starboard  vortices  are  labeled  by  the  letters  P and  S respec- 
tively while  the  symbol  (*)  refers  to  a single  vortex  whose  origin  is  undefined. 
In  Figs.  5-1  through  5-10  the  LDV  system  consistently  locates  the  port  and 
starboard  vortices.  At  late  times,  when  lateral  drift  of  one  vortex  out  jf  the 
LDV  field  of  view  can  occur,  discrimination  of  the  port  and  starboard  vortices 
is  limited. 

The  lateral  trajectory  of  the  wake  vortices  is  indicated  by  the  height 
versus  y position  and  vortex  age  versus  y position  plots.  For  the  sample 
runs  shown  in  Figs.  5-1  through  5-10  the  wake  vortex  remains  in  the  region 
+150  ft  about  the  runway  centerline  for  20  to  60  sec.  Some  scatter  can  be 
seen  in  the  lateral  trajectories  attributed  to:  (1)  the  actual  vortex  wandering; 

(2)  the  range  resolution  of  the  LDV  system;  and  (3)  the  resolution  of  the  vortex 
track  algorithm.  Filtering  the  data  or  a least  squares  fit  to  the  vortex  tra- 
jectory could  eliminate  some  of  the  scatter.  However,  the  existence  of  the 
wake  vortices  in  the  vicinity  of  the  runway  can  be  clearly  determined  from 
the  measurements. 

In  addition  to  the  scatter  in  the  vortex  track  measurements  attributed 
to  ambient  conditions  and  to  computer  processing  techniques,  some  of  the 
scatter  resulted  from  constraints  in  the  LDV  system  hardware.  Two  possible 
sources  of  error  in  the  LDV  system  included:  (1)  error  in  vortex  track  loca- 
tion due  to  actual  versus  commanded  focus  position,  and  (2)  failure  to  record 
vortex  signal  during  periods  when  the  onboard  computer  was  dumping  the  filled 
up  buffer.  Since  the  location  of  the  LDV  focal  volume  was  determined  from  the 
commanded  rather  than  the  observed  mirror  orientation  and  since  a small  lag 
occurred  between  the  commanded  and  observed  mirror  position,  the  measured 
vortex  location  could  be  in  error  by  +20  ft  for  the  typical  finger  scan  rates 
selected  at  JFK.  Hardware  problems  with  the  computer  interrupt  dictated  a 
single  buffer  data  collection  mode  on  the  onboard  computer  which  resulted  in 
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^ data  being  rejected  whenever  the  computer  was  in  the  process  of  dumping 

the  filled  up  buffer.  As  a result,  intermittent  skips  in  the  recorded  vortex 
tracks  lasting  a fraction  of  a second  could  be  observed  during  the  measure- 
ments. However,  it  is  noted  that  the  above  hardware  problems  did  not 
i seriously  hamper  the  capability  of  the  LDV  to  identify  and  to  monitor  wake 

vortices  in  the  vicinity  of  the  runway.  Appropriate  modifications  can  be 

made  to  the  LDV  system  to  circumvent  the  cited  hardware  limitations. 

f 

i 

; The  vertical  trajectory  of  the  wake  vortices  is  indicated  by  the  vortex 

^ age  versus  height  and  height  versus  y-position  plots.  For  the  sample  runs 

^ shown  In  Figs.  5-1  through  5-10  the  vortex  pair  descends  at  approximately 

constant  velocity  until  it  is  within  50  to  100  ft  of  the  ground.  The  scatter  in 
the  vortex  altitude  measurements  appears  to  be  larger  than  in  the  vortex 
lateral  position  measurements.  This  is  attributed  to:  (1)  the  wake  vortices 
can  become  asymmetric  during  the  descent  process,  and  (2)  the  uncertainty 
in  the  altitude  is  related  to  the  product  of  the  uncertainty  in  the  elevation 
angle  and  range. 

Decay  of  the  wake  vortex  peak  velocity  is  illustrated  by  the  vortex  age 
versus  peak  line -of-sight  velocity  plots  in  Figs.  5-1  through  5-10.  The  peak 
• line-of-sight  velocity  is  a measure  of  the  maximum  rotational  velocity  of  each 

j vortex.  In  general,  no  significant  decrease  of  the  vortex  rotational  velocity 

t occurs  initially  in  the  vortex  wake,  i.e.,  0 to  40  sec  after  aircraft  passage, 

while  at  later  times  the  peak  velocity  decreases  markedly.  Some  scatter 
' can  be  noted  in  the  peak  velocity  time  histories  which  is  attributed  to:  (1)  the 

different  velocity  decay  rates  for  the  port  and  starboard  vortices,  (2)  unsteadi- 
ness in  the  vortex  flow  and  the  lack  of  axial  symmetry,  and  (3)  the  uncertainty 
In  determining  the  exact  vortex  centroids  due  to  spatial  resolution  and  lack  of 
sample  points.  However,  the  velocity  decay  trends  shown  in  Figs.  5-1  through 
5-10  are  in  agreement  with  the  plateau  and  the  1/t  or  i/VT  type  of  vortec  decay 
characteristics  observed  in  wind  tunnel  and  water  tank  studies  by  others  (Ref.  7). 

The  wake  vortex  trajectories  measured  with  the  LDV  system  at  JFK  are 
given  in  Appendix  A and  show  the  same  trends  as  the  sample  plots  discussed 


6.  RESULTS  OF  WIND  MEASUREMENTS 


The  three-dimensional  wind  field  was  monitored  as  a function  of  altitude 
near  the  middle  marker  position  of  runway  31R  at  Kennedy  International  Air- 
port several  times  daily  during  the  LDV  field  tests.  The  objective  of  the  wind 
data  collection  was  to:  (1)  demonstrate  the  capability  of  the  Lockheed-Huntsville 
LDV  for  measuring  ambient  atmospheric  wind  fields;  (2)  compare  the  LDV 
measurements  with  wind  measurements  obtained  by  conventional  anemometers 
at  the  test  site;  and  (3)  from  the  above  data,  establish  the  operational  capa- 
bilities, resolution,  and  integrity  of  the  LDV  for  wind  monitoring  at  terminal 
areas. 

During  the  JFK  field  tests,  a considerable  amount  of  wind  data  was 
collected  with  the  LDV  as  shown  in  the  data  logs  in  Appendix  A.  However, 
a detailed  comparison  of  the  LDV  wind  measurements  with  on-site  wind  meas- 
urements from  the  DOT-TSC  meteorological  towers  was  hampered  by  lack 
of  data  from  the  meteorological  towers.  In  subsequent  tests,  summarized 
in  a separate  report,  LDV  u,  v,w  measurements  were  carefully  correlated 
with  instrumented  tower  measurements;  the  results  showed  that  the  LDV  is 
a highly  sensitive  wind  monitoring  device  (Ref.  6). 

Wind  profile  measurements  were  carried  out  at  JFK  with  the  LDV  sys- 
tem operating  in  the  VAD  mode.  The  processed  wind  speed,  direction,  and 
crosswind  and  downwind  profile  are  shown  in  Figs. 6-1  and  6-2  for  sample 
cases.  For  comparison,  the  hourly  20-foot  level  NOAA  groundwind  meas- 
urements at  JFK  are  also  plotted  for  each  of  the  wind  profiles  as  obtained 
from  the  surface  weather  charts.  Good  correlation  is  noted  between  the 
Lockheed-Huntsville  LDV  measurements  fitted  by  a least  squares  curve  fit 
and  the  ground  wind  observations  in  Figs.  6-1  and  6-2.  A power  law  wind 
speed  proftle  and  a linear  wind  direction  profile  were  fitted  to  the  LDV  VAD 


measurements  at  different  altitudes  so  as  to  minimize  the  sum  of  the  squares 
of  the  derivations  of  the  given  points  from  the  curve.  The  LDV  measurements 
have  also  been  compared  with  observations  from  the  DOT-TSC  instrumented 
meteorological  tower  at  JFK.  The  LDV  system  was  located  3,700  and  850  ft, 
from  the  DOT  135  and  40  ft  towers,  respectively.  The  NOAA  tower  was  located 
approximately  2 miles  from  the  LDV  system.  The  results  in  Figs.  6-3  and  6-4 
show  the  wind  speed  measurements  obtained  from  the  LDV  VAD  scans  and  the 
meteorological  towers.  The  measurements  in  Fig.  6-5  show  the  line-of-sight 
velocity  obtained  with  the  LDV  system  focused  at  the  location  of  the  anemometer 
on  the  40'  tower.  The  wind  measurements  from  the  40'  DOT  tower  and  the 
20'  NOAA  tower,  resolved  about  the  LDV  line-of-sight,  are  shown  in  Fig.  6-5 
for  comparison.  Again,  good  agreement  is  noted  between  the  LDV  and  meteoro- 
logical tower  measurements. 


Wind  ProfUe,  JFK  Runway  31R  (10/30/75) 


Runway  Wind  Velocity  (ft/sec) 


Fig.  6-2  - Wind  Profile,  JFK  Runway  31R  (10/30/75) 
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Fig. 6-3  - Comparleon  of  Meteorological  Tower  Wind  Speed  Meaauremente 
with  Lockheed -Huntsville  LDV  at  JFK  on  11/5/75  14:06-14:36 
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Comparison  of  LDV  Line  of  Sight  and  Corresponding  Meteorological 
Tower  Wind  Speed  Measurements 


7.  CONCLUSIONS 

The  operational  capability  of  the  Lockheed-Huntsville  mobile  laser 
Doppler  veiocimeter  system  for  the  measurement  of  winds  and  wake  vortices 
at  terminal  areas  has  been  demonstrated.  Wake  vortex  measurements  have 
been  obtained  under  a range  of  operating  conditions  and  the  existence  of  wake 
vortices  in  the  approach  corridor  was  monitored  successfully  during  the  field 
tests.  A useful  data  base  of  wake  vortex  trajectories  were  collected. 

Analysis  of  the  performance  of  the  LDV  suggested  that  the  following 
modifications  could  further  improve  the  remote  sensing  capabilities  of  the 
system: 

1.  Incorporate  actual  rather  than  commanded  elevation  angle  signal 
into  the  scanner  and  double  buffer  the  onboard  computer  to  elimi- 
nate anomalies  in  the  data  acquisition. 

1.  Increase  the  data  acquisition  rate  and  resolution  of  the  system  by 
integrating  a filter  bank  into  the  signal  processor.  — 

3.  Explore  different  scan  configurations  and  different  vortex  dis- 
crimination concepts  to  improve  the  definition  of  the  vortex  tracks. 

The  results  of  the  research  program  obtained  a valuable  data  of  wake 
vortex  trajectories,  demonstrated  the  basic  reliability  of  the  LDV  system, 
and  suggested  techniques  for  refining  the  capability  of  the  LDV  remote 
sensing  system. 
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